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ABSTRACT 


Methods  are  explained  for  deriving  the-  differential  equations  des¬ 
cribing  a  general  .missile  airframe  flying  In  three-dimensional  space. 

All  the  forces  and  moments  acting  on  the  missile  are  written  as  vectors. 
These  vectors  axe  readily  expressed  In  suitable  coordinate  systems. 

After  they  Axe  so  expressed,  they  are  transformed  to  a  common  coordinate 
system  where  the  solutions  are  carried  out.  The  desired  output  informa¬ 
tion  Is  then  transformed  to  a  fixed  system  an  the  ground  and  presented  in 
texms  of  that- system.  A  short  exposition  of  vector  analysis,  to  the  ex¬ 
tent  required.  Is  included,  -sad  matrix  notation  is  emphasized  throughout. 

Methods  (utilizing  vector  analysis  end  coordinate  system  transforma¬ 
tions)  are  presented  for  handling  thrust  misalignment,  the  general  moment 
of  inertia  tensor,  curvature  of  the  earth,  rotation  of  the  earth,  targets, 
etc. 

The  methods,  as  t ^ployed;  generate  an  exact  representation  of  the 
flying  missile;  however,  wherever  possible,  approximations  are  suggested 
which  simplify  the  equations.  Errors,  not  only  from  approximations,  but 
from  other  sources,  are  pointed  out  when  such  situations  arise. 

Since  there  are  so  many  types  of  guidance  philosophies,  only  very 
general  statements  are  made  concerning  guidance.  However,  the  guidance 
philosophies  for  two  speclfio  types  of  missiles  are  included  as  examples. 
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FvPMMEfEAIS  FCrt  THE  GENERATION  OF  A  MATH  MODEL  FOR  A  MISSILE  INTHREE- 

EOMESSICKAL  SPACE  WITS  3 IX.  DEGREES  OP  FREEDOM. 

pmOKJCTION 


The  following  series  of  notes  give 3  the  derivation  for  the  general 
equations  'which  describe  any  missile  system  in  three  dimensional  space. 

Many  people  have  expressed  a  desire  to  have  these  fundamentals  expounded.  In 
detail  and  also  have  the  derivations,  themselves,  in  printed  form  for  future  - 
reference  hence  this  series. 

These  notes,  which  will  al3o  be  given  as  lectures,  have  previously 
been  given  to  a  select  group  in  a  sort  of  "dry  run".  That  group  fait  that 
these  notes  would  provide  a  valuable  reference,  and  could  be  used  also  as  a 
tiindaxental  basis  In  deriving  equations  for  any  missile  system,  since  the 
equations  for  a  general  missile  system  are  of  the  same  type  no  matter  what 
type  of  missile  is  used.  The  airframe  equations  (which  these  lectures  are 
primarily  concerned  with)  are  more  or  less  fundamentally  identical  In  every 
csss  of  a  flying  missile.  We  are  dealing  with  a  self -powered  free-free  beam  ■  ' 
in  space  which  is  acted  upon  by  aerodynamic,  propulsive,  and  gravitational 
forces.  No  matter  what  type  of  missile  ve  are  using  the  same  forces  will  be 
acting  on  it.  The  method  for  obtaining  the  accelerations  and  trajectory  follow 
the  basic  physics  of  the  free-free  beam  configuration. 

Of  course,  by  three  dimensional  space  we  mean  the  atmosphere  which  ‘ 
surrounds  us.  The  phrase  "six  degrees  of  freedom"  describes  the  number  of 
equations  required  to  solve  this  system:  three  for  the  solution  of  the  for¬ 
ces  and  three  more  which  give  solutions  for  the  moments.  From  these,  Ve  axe 
able  to  obtain  enough  independent  equations  to  determine  completely  the  un¬ 
known  quantities  concerning  the  missile.  The  term  free-free  beam  mentioned 
earlier  simply  describes  the  appearance  of  a  flying  missile. 

We  will  begin  with  the  very  basic  fundamentals  and,  stop  by  step,  accrue 
the  necessary  terms  to  completely  describe  the  missile  system.  Ba  many  cases 
the  Information  which.  Is  contained  in  this  document  may  be  something  which  the 
afier  knows  veil  and  it  will  be  a  redundancy.  Hrrarer,  for  many  others  who 
have  not  had  ths  opportunity  to  meet  this  exposition,  this  will  not  be  the 
caue. 

The  basic  organisation  will  proceed  from  the  development  of  the  force 
and  moment  equations  and  the  development  of  the  coordinate  systems  necessary 
to  describe  the  forces  and  moment q,  to  the  matrix  transformation*  which  will  he 
employed.  Along  the  vay  ve  will  introduce  the  concept  of  relative  wind  end 
methods  by  which  ve  may  study  such  things  as  thrust  misalignment.  The  foogd 
lecture  period  will  he  devoted  to  describing  the  aspects  of  vector  analysis 
which  ve  will  seed  for  our  problem*.  If  the  reader  or  listener  is  familiar, 
with  this  vector  algebra  he  may  skip  this  part  with  no  lots  of  understaudfiBg**  ‘ 
It  Is  Included  here  primarily  to  review  and  remind  us  what  part  of  the  large 
body  of  knowledge  called  vector  analysis  ve  will  use. 
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Whenever  ve  with  to  solve  equations  concerning  any  missile  systes  ve 
ere  primarily  Interested  In  determining  bow  high  or  hew  far  tbs  missile  will 
go*  Tbe  final  result  ve  wish  to  obtain  from  any  series  of  equations  are  the 
ranges  and/or  velocities  at  any  Instant  of  tine*  If  we  can  obtain  expressions 
which  detexalne  tbe  position  of  tbe  missile  with  respect  to  soon  fixed  place, 
we  will  bare  tbe  solution  to  alnost  all  prbblsns  with  which  ve  will  be  dealing; 
certainly,  those  concerned  with  trajectory  Information* 

VS  will  begin  by  describing  tbs  forces  acting  on  tbe  nlsslle*  To  do  this, 
ve  will  solve  Hevton1*  so-called  second  lav  of  notion,  tbe  fore?  equation.  In 
vector  notation  this  Is: 

x 

where  tbe  vector  £F  is  tbe  total  force,  n  Is  tbs  ns ss  of  tbs  nlsslle  at  any 
instant,  end  vector  a  Is  tbs  eooaleratlon*  Vs  can  lanadlattly  break  this 
eqostlen  dons  Into  tbs  following  form. 

£fz  py  f  Fa  ¥fr  + 

In  thdLs  expression  Fjr  Is  tbe  gravitational  force,  tbe  next  tern  Is  the  aaro- 
dynanlc  force,  fonowed  by  tbe  thrust  fores.  So  neks  things  complete,  a  tern 
Is  put  In  called  F#t«m  In  oase  our  nlsslle  has  soon  odd  type  of  force  which 
is  not  contained  In  tbe  gravitational,  aerodynamic,  or  thrust  terns.  For  ex¬ 
ample,  •  rotating  nlsslle  experiences  s  side  force  celled  tbe  Megans  Kffeot— 
a  force  Which  would  be  izduded  In  .  In  general  tbe  only  forces 

acting  on  tbe  nlsslle  ere  the  grevitetloral  fores  e,  the  aerodynamic  force,  end 
tbe  power  or  thrust  force.  Throughout  these  lectures  ve  will  issue*  that  there 
It  no  "other''  force.  Such  forces  as  Coriolis  force  and  centrifugal  farce  ere 
fictional  forces  which  nsrely  describe  the  transformation  between  e  non -accel¬ 
erated  coordinate  system  and  one  which  is  accelerated.  In  this  derivation 
there  will  be  no  necessity  for  ever  exsairlrg  such  "fictional"  forces  spedflv- 
oally.  If  vs  wish  to  actually  look  at  then,  they  can  be  found  In  tbe  equations, 
but  they  arise  naturally  out  of  tbs  transformations  which  we  will  use  In  de¬ 
riving  tbe  equations. 

Vevton's  nsc  end  lev,  as  stated  above,  Is  only  true  (disregarding  ralatlv* 
Istio  effects)  In  an  “inertial"  system.  By  an  "Inertial"  system,  ve  mean  one 
Which  Is  not  accelerated.  Therefore,  we  must  designate  sons  plaoe,  preferably 
where  we  are  steading,  as  being  unaooal stated.  Oeaerally,  this  inertial  sys¬ 
tem  is  put  on  tbe  grot end,  sad  so  la  this  series  we  are  going  to  anil  our 
insrtlal  system  tbs  "groual  system*.  Actually,  for  a  oomplsta  and  aocurate 
dssorlptlon  of  tbe  flying  mlsslls,  we  should  put  our  Inertial  systea  at  tbs 
a  art  Mr  of  tbs  earth,  because  tbs  earth  rotates  cm  Its  sols* 


In  this  csss,  the  earth  would  rotate  around  the  "inertial"  system  which  re* 
jat-ins  fixed  at  the  center  of  the  earth.  The  point  of  launch,  than,  would 
rotate  with  the  earth  with  velocity  equal  to  the  tangential  velocity  of  the 
earth's  surface.  However,  for  ease  of  derivation,  we  are  going  to  assume 
that  the  earth  does  not  rotate.  (Since  the  rotation  Is  slow,  in  many  cases 
the  induced  error  is  small.  enough  to  be  ignored.)  If  it  becomes  necessary 
(and  it  would  if  we  were  working  with  long  range  missiles)  ws  could  put  our 
Inertial  system  at  tbs  canter  of  the  earth  and  the  transformations  to  he 
subsequently  described  would  be  of  exactly  the  ease  type.  SirCtthls  series 
of  lectures  is  intended  to  Inpart  Information,  ws  will  not  go  Into  the  da* 
tail  of  the  transformations  frost  the  center  of  the  earth.  She  only  thing 
that  this  type  of  transformation  would  do  would  be  to  give  more  complexity 
to  the  resulting  equations  but  add  nothing  to  the  theory.  Another  thin; 
that  va  are  going  to  assume  is  that  the  earth  la  flat.  Generally,  the 
error,  because  of  this,  le  not  greet  if  the  trajectory  is  relatively  short. 

If  the  trajectory  ware  long,  it  would  require  that  va  make  corrections  to 
oux  equations  to  take  into  account  the  curvature  of  the  earth.  Again,  this 
merely  makes  the  resulting  equations  more  complicated  and,  since  it  adds 
nothing  particularly  to  the  theory,  wa  will  not  go  into  title  aspect  of  the 
problem. 

He  will  begin  by  issued  ng  our  inertial  coordinate  system  to  he  fixed 
to  the  earth  at  the  plaoe  of  launch  of  the  ml  sal,  it.  She  inertial  system,  or 
the  ground  system  as  we  will  oall  it,  u  a  right -tended  orthogonal  aMteeiea 
eat  of  ana.  They  will  be  defined  by  three  unit  vectors  jr,  and X, 
Before  we  can  proceed  any  further  uo  must  define  what  we  mean  By  positive  and 
negative  translations  end  rotations.  If  wa  translate  in  a  positive  direction, 
we  move  from  the  oenter  or  origin  of  the  coordinate  syvten  outwards  la  the 
direction  which  a  nit  vector  points.  She  follomlag  ihgnm  will  lUntsmte 
what  wa  moan  by  positive  rotation} 


If  ve  rotate  about  a ay  of  tie  three  unit  vector*,  or  exes,  then  the  rotation 
nill  be  positive  if,  when  va  are  rotating  about  the  first  unit  vector,  g^,  the 
rotation  is  from  unit  vector  number  two  toward  unit  vector  number  three*  If 
we  rotate  about  the  second  unit  vector,  the  rotation  is  positive  if  rotates 
toward  "g,,  Rotation  la  positive  about  the  third  vector  if  the  rotation  is  in 
the  tense  from  gjj.  to  g£.  therefore,  whenever  we  rotate  if  ve  always  rotate  in 
these  directions,  that  is,  from  ana  to  two,  two  to  three,  and  three  to  one, 
then  our  rotations  will  be  positive,  and  all  tbs  agnations  will  automatically 
have  the  correct  sign. 

We  are  now  In  e  position  to  further  define  'the  forces  in  the  previous 
equation*,  far  sample  the  gravitational  fores  always  points  to  the  center  of 
the  earth*  If  we  assume  that  our  third  ground  axis,  that  is,  "83,  is  pointing 
upward  from  tbs  center  of  the  earth  then  the  gravitational  force  is  given  as: 

Fjr  5  “*7,133s 


The  aerodynamic  forces,  on  the  other  hand,  arc  related  specifically  to  the 
direction  which  the  aiaai.le  is  flying*  It  would  be  extremely  difficult  to 
describe  the  aerodynamic  farces  in  terms  of  the  ground  system  which  we  have 
just  defined  above*  Therefore*  for  this  purpose,  ve  will  define  a  new  coord¬ 
inate  system  celled  the  wind  system*  It  is  another  orthogonal  cartesian 
coordinate  system  with  three  unit  vectors  %  and  v,.  In  this  system  ve 
writs  down,  exactly,  the  aerodynamic  forces  in  our  equations*  These  are 
written  set  ^  ^ 

ft*  D&1  + 


m  the  wind  eyetem  the  terms  used  for  the  aerodynamic  force  ere  dreg,  aide 
force,  and  lift,  and  these  tarns  refer  only  to  sarodynamlo  forces  defined  in 
the  wind  system.  This  wind  system  is  so  set  up,  that  the  first  unit  vector 
%  always  points  in  the  direction  of  the  relative  wind  vector.  The  other 
two  vectors  ^  and  arm  such  that  ve  can  rotate  through  two  angles,  called 
angles  of  sttaak*  ana  these  veatoro  will  be  oriented  with  the  body  of  the 


angles  of  attack, 
missile. 


The  thrust  force,  IL,  will  be  defined  in  a  system  called  tbs  thrust  jys- 
tsm*  The  coordinates  will  be  designated  by  three  orthogonal  unit  rectors  p^ 
"$2,  and  P3,  ve  choose  "p"  rather  than  "t"  for  thrust,  beosuse  *t"  will  to*  used 
as  temperature,  tine,  and  reference  to  the  targst  in  other  equations.  Actually 
we  can  oaU  V  tbs  'push*  of  thm  miseil*.  We  will  dmearibe  the  thrust  system, 


that  is,  thm  R »*,  so  that  the  forum  doe  to  thf  ejection  of  material  from  the 
ml  mails,  the  thrust,  acts  in  thm  direct! on  of  p*.  Than  ve  can  simply  rotate 
through  two  angles  in  order  tc  Una  up  our  thrust  system  with  thm  admails  body* 
If  we  asmtsM  there  is  no  thrust  misalignment,  than  thm  angles  between  the  "p" 
system  and  tbs  mdaatle  body  will  toe  saro  —  and  the  thrust  will  toe  pr  aetmed  to 
act  through  the  miemile  osmtearllne.  If  we  ever  wish  to  investigate  thrust  mlm- 
aligamaat,  ve  will  already  have  the  math— ties  for  doing  so. 


Ve  will  mead  to  ftaasritoa  one  other  coordinate  system  whose  ajgrdfioanoe 
mill  toaocma  more  apparumt  later.  _  Ode  is  called  the  body  system,  end  is  desig¬ 
nated  by  the  three  wait  vmotorm  JL,  amd>».  Vs  arm  actually  going  to  solve 
all  of  osar  eqaatVyom  in  this  body  system.  Bos  system  is  defined  in  the  follow¬ 
ing  nsarnan  the  ev,  that  is  ths  body-1. axis  is  aligned  so  that  it  is  slung  tbs 
oenter  Una  of  the  aissils*  Ha  wfU  assume  that  the  origin  of  this  system  is 


at  the  center  of  gravity  of  the  missile.  The  other  two  vectors,  bo  and  bj,  can 
be  defined  more  or  less  at  the  analyst's  discretion*  However,  it  is  apparent, 
or  at  least  will  be  apparent,  that  certain  types  of  orientations  are  better  than 
others. 

The  following  diagram  illustrates  the  coordinate  systems  which  have  Just 
been  described:  those  coordinate  systems  which  we  nailed  the  ground  system,  the 
body  system,  the  thrust  system  and  the  wind  system* 


VECTORS 


Unit  Vectors 


In  the  work  vhich  we  are  concerned  with,  that  la  the  description  of 
missiles,  their  trajectories,  and  responses,  the  most  powerful  mathematical 
tool  ve  car  use  is  rector  analysis.  In  any  specs,  say  two-dimensional  or 
three-dimensional  in  vhich  ve  are  generally  concerned  in  this  type  of  study, 
vo  can  define  the  position  of  any  point  by  means  of  two  lines  on  a  plane  or 
by  three  lines  in  a  three-dimensional  space,  so  long  as  those  lines  do  not 
coincide.  For  example,  if  ve  draw  tvo  lines  in  a  plane  as  in  the  following 
diagram^ 


then  we  can  get  to  any  point  on  this  plana  by  moving  only  parallel  to 
these  two  lines.  In  the  above  diagram  we  vest  negatively  along  A  and  than 
vent  positively  along  B,  where  positive  is  in  the  direction  of  tbs  arrows* 

If  ve  aet  up  a  suitable  scale  on  the  lisas  A  and  B  so  that  they  hove  a 
definite  length,  and  if  we  decide  that  they  point  in  a  certain  direction, 
then  these  are  called  vectors*  If  the  lir.cn  A  and  B  were  of  unit  length  (in 
any  measuring  system)  they  are  oalled  unit  vectors.  Furthermore  if  ve  made 
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these  vectors  orthogonal,  that  Is  at  right  angles,  as  la  the  following  diagram: 

f 

3 - *8 


then  these  two  vectors  of  length  one  and  at  right  angles  to  each  other  are 
called  orthcooxmL  vectors-  In  three  dimensional  space  we  will  draw  one  acre 
vector  of  wit  length  which  Is  orthogonal  to  both  1  and  ?  and  we  will  have  a 
three  dimensional  system  with  the  third  vector  &  We  can  now  describe  any 
point  In  space  with  respect  to  these  unit  vectors*  Thus  the  vector  ?  in  the 
following  diagram  would  be  written  In  this  system  as:  P  «=  x3f  +  yB.4  zC, 


which  simply  tells  as  to  travel  x  units  along  X  y  units  along  "B,  and  %  units 
along  C  to  arrive  at  point  P. 

^  'lolng  back  to  our  original  ground  system,  we  have  three  unit  vectors, 
X,  g £,  and  g’j,  mutually  orthogonal: 


These  vectors  form  an  orthononal  coordinate  system*  We  are  also  going  to 
restrict  ourselves  to  right-handed  coordinate  systems*  A  right-handed  coord¬ 
inate  system  la  defined  as  follows: 

The  direction  of  Is  the  direction  of  advance  of  a  right-handed  screw 
when  It  Is  rotfted'Trom  %  towards  %•  The  direction  of  Hk  Is  the 
direction  of  advance  of  aright -banded  screw  when  rotated  from  to¬ 
wards  .  The  direction  of  ^  Is  the  direction  of  advance  of  a  right- 
handed  screw  when  rotated  from  gj  towards 

It  Is  possible  to  use  a  left-handed  coordinate  system.  Indeed,  In  Russia  the 
loft-handed  coordinate  system  Is  tbs  nun.  However  throughout  the  United  States, 
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England  and  much  of  Europe  the  standard  la  the  rigirs  “banded  coordinate  spates, 
it  is  the  one  ve  vill  use  hare,  and  it  is  the  one  in  ifeleh  it  is  most  easy  to 
visualize  our  functions. 

Assume  ve  have  defined  an  orthoscaaaal  right -handed  coordinate  system  and 
ve  have  scae  vector  in  space  called  vector  X* 

t 


1 - 

In  the  above  diagram  X  has  sene  length  and  it  is  oriented  in  a  certain  way 
with  respect  to  the  coordinate  system#  7c  simplify  our  analysis  ve  vill 
only  talk  about  two  unit  vectors  gx  and  "ga*  If  the  length  of  vector X  is 
J]A"  units  and  tbs  angle  between  a  line  parallel  to  the  gL  axis  and  the  vector 
A  is  0,  tinsn  tho  length  of  the  line  parallel  to  the  Xx  aJns  is  e^ual  tolXjooe 
9,  and  the  length  of  the  lige  parallel  to  the  axis  isiXjsin  8^  Therefore, 
ve  can  aay  t^at  the  vector  A  is  j  where  the  abso¬ 

lute  value,  tkl  represents  the  magnitjDie  of  the  vector  without  any  indication  of 
direction.  The  direction  of  vector  X  is  given  by  Xx  and  g^. 

Addition  and  Subtraction  of  Vectors 

Assume  ve  have  two  vectors  T  and  5  which  have  components  la  the 
system  as: 

A-  At,  +At% 

S=  3\%+BtSt  *5jJj 

These  vectors  con  be  manipulated  by  means  of^veotor  algebra.  For  instance, 
two  vectors  can  be  added.  The  addition  of  A  and  "a  is  defined  as  follows: 

A+B  r  (Afr&lSS'  +  (fi*+8s$%  + 

The  zjesult  of  adding  the  two  vectors  A  and  B  is  that  ve  add  the  component*  of 
the  vector*  together  (the  distance  ve  go  in  the  direction  of  XlO  Similarly, 
ve  add  the  components  of  end's**  Ve  can  demonstrate  the  addition^  thage 

two  vectors  as  in  the  following  diagram  (assuming  T  and  "I  only  have  and  g$ 
cuapcmants) : 
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Zt  can  te  easily  seen  that  adding  the  components  of  the  vectors  A  and  B,  to 
each  other  is  tantamount  t'o  adding  the  two  vectors  with  their  respective  direc¬ 
tions  tail  to  head  and  drawing  a  line  between  the  tail  of  the  first  to  the  head 
of  the  second.  She  above  diagram  illustrates  another  important  principles  the 
actual  position  in  space  of  those  vectors  is  immaterial.  Ve  can  move  then  any¬ 
where  we  want  as  long  as  we  keep  their  directions  and  magnitudes  the  same,  Tfcus 
we  could  move  the  B  vector  up  and  attach  it  to  the  A  vector  to  produce  the  re¬ 
sultant  vector.  The  resultant  of  adding  these  tw  vectors  is  a  third  vector, 
jtf  B.  Subtraction  of  vectors  would  be  accomplished  simply  by  changing  the  sign 
of  the  components!  *  As  can  be  seen  by  a  simple  sketch,  this  is  the  same  as  if  ve 
reversed  the  direction  of  the  subtracted  vector  sad  did  not  interfere  with  the 
magnitude: 


The  Pot  Product: 

Mcltiplioctica  of  vectors  can  be  perfoxmed  in  two  different  ways.  One 
of  these  is  the  dot  product  or  scalar  product,  the  other  is  the  cross  or  vec¬ 
tor  product.  The  dot  product  gives  a  scalar  -.umber,  that  is,  a  number  having 
a  magnitude  but  no  direction.  The  cross  product  gives  a  new  vector.  In  each 
case  we  perform  a  type  of  multiplication  between  two  vectors  or  perhaps  among 
several  vectors.  The  dot  product  can  be  defined  in  two  ways.  Using  the  vec¬ 
tors  "X  and  it  described  above  ve  can  define  the  dot  product  between  them  as: 


where  9  is  the  angle  between  the  vectors  A  and  B.  If  ve  write  the  vectors  A 
and  a  in  their  component  fora,  then  the  dot  product  between  the  two  is  the 
regular  seals?  or  number  product  of  each  component  and  then  the  products  of 
these  components  are  dimmed  together  as  balov: 


B*S  B\ 

A*lSz/lt3,+Aj8i3 


Xu  the  type  of  analysis  which  we  plan  on  doing,  the  most  Important  type 
of  dot  product  is  that  between  unit  vectors*  going  back  to  the  first  definition 
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let  us  take  a  dot  proiuot  as  va  would  a  regular  algebraic  product  between  two 
polynomials .  Thus : 


At  this  point  ve  should  mention  that  the  scalars  Multiplying  rectors  can  be 
factored  and  moved  to  the  front  of  the  expression.  Thus  the  first  term  in 
the  above  equation  can  be  written  as  *\LB,  g,  ."g  .  •  We  know  ignr  to  Multiply  real 
numbers,  bo  vc  Must  find  the  dot  product  between  and  gj_.  In  each  of  the 
above  terms  if  vw  use  the  first  definition  of  the  dot  product  we  see  that  we 
have  the  following  form  for  the  product  of  these  two  vectors: 

We  see  that  the  dot  product  between  %  and  product  of  the  magnitudes 

of  g-_  times  the  cosine  of  the  angle  between  them.  In  this  case  the  magnitudes 
are  coe  times  one  and  the  angle  between  a±  and  gj.  is,  of  courser  aero.  There¬ 
fore,  “g,  ."g^  is  equal  to  one.  On  the  other  hand  the  product  g±*&\  is  com  times 
cue,  times  the  cosine  of  the  angle  between  g|  and  *&•  We  have  defined  our 
coordinate  system  so  the  %  is  orthogonal  to  gu  and  tberefocre,  the  angle  be¬ 
tween  g^  and *g^  is  90°.  She  cosine  of  90°  is  zero.  Thus,  eod 

equal  zero.  The  only  terms  loft  are  the  first  and  last  terms  and  the'ViXue  of 
the  dot  product  Is  AjBvf  AoB.,.  This  allows  us  to  state  In  words  a  rule  for  the 
dot  products  of  vectors  when  defined  In  orthonoraal  systems:  "The  dot  product 
between  unit  vectors  equals  one  if  the  unit  vectors  are  the  same  and  It  equals 

zero  if  they  are  not  the  same." 

* 

The  Cross.  Product 

The  definition  of  the  cross  product  is: 

A  X&  =  * 

That  l^the  cross  produat^of ~A  end’s  Is  the  scalar  product  of  the  Magnitude  of 
"X  times  the  magnitude  of  B  times  the  sine  of  the  angle  between  thsm.  The 
direction  £  is  a  unit  vector  orthogonal  to  the  plane  of  Tf  and  B.  Thus: 

1ix£~0fc,net 


1 r 


Furthermore,  the  direction  of"?  is  pointing  in  the  direction  of  the  advance  of 
a  right-handed  screw  if  we  rotate  X  towards  X. 

*  The  fact  that  we  have  a  vectorial  result  is  not  an  accident,  but  the  proof 
relies  on  tensor  analysis  and  will  not  be  dealt  with  here. 
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B*  protect  of  B  cross  A,  then,  would  bo  la  tbs  direction  of  tbs  advance 
of  e  right-handed  screw  If  ve  rotate  B  towards  X*  Bis  can  be  seen  tram  the 


Bus  vs  can  say  that: 

See  arose  product  between  unit  rectors,  which  ve  will  generally  be  concerned 
with;  is  as  follows: 

iron  tbs  definition,  :,  If  A«8  tben<9e0  m&’tj&S* 

Bw  cross  products  of  the  ^  vectors  are:  i 


31  *},*%*& *t**%s° 


As  ve  atated  above,  tbs  cross  product  produces  a  new  vector  pointing  In  the 
direction  of  tbs  advance  of  a  right-handed  screw,  If  ve  rotate  the  first 
vector,  la  the  teas,  toward  tbs  seoond*  Tram  tbs  definition  at  positive 
w rotation  and  orthonoamJlty,  the  above  definitions  for  the ^,'s  are  laned- 
late* 


A  convenient  aetbod  of  deteralning  the  cross  product  of  two  veotors 
is  by  using  the  cetcndnaat  defined  below: 

£<Va*  $ m +/)&+*/£  j  3e*'$ f 
ThM  |  %,%zO\tnrA> 

0t  Bt  & 


xo 


The  two  definitions  for  the  cross  product;  that  is,  the  determinant  . 
method  and  the  de  flni  t  i  01^  /fy  £  zf/fj-f  ffjsi  n  &  S»  are  completely  equivalent* 

The  proof  of  this  is  easily  shown  by  expansion  and  will  not  be  given  here. 

Throughout  our  wori  in  the  analysis  of  missile  systems  we  will  be 
using  the  dot  product  bet-seen  unit  vectors  and  tbs  cross  product  between 
unit  vectors.  This  is  because  we  will  always  describe  our  forces,  our  mo¬ 
ments  and  other  aspects  of  the  response  of  the  missile  is  terms  of  magnitude, 
with  the  directions  indicated  by  the  unit  vectors. 

Triple  Cross  Product  and  triple  Scalar  Product 

From  time  to  time  vector  products  of  the  form  '&(&£)  arise.  Such  a 
vector  is  called  the  triple  cross  product  or  triple  vector  product.  That  tha 
triple  cross  product  is  hot  necessarily  associative  is  very  important.  That 

This  is  very  easy  to  show.  Suppose  that 

S’-  8i$,+3£§2+/bijb 
$*Clt+-Cjtf,+C3  % 

Than  to  solve  the  product  !X(M);  we  first  find: 

txc-  l  $2  Si  -  -B^tt 

8.  82 

C>i  Cg  c$ 

k>KieJ*  Jivcss  % 

TSx$x'0=  a  fc  * 

BgCj-EbCg  BjCr&tCa  BiCyfaCi 

Solving: 

~Sx  (3*C)  =  [Si  (m  StC,)  -Bj  (BA  '&  $ +[<b(etC,-bcd-6lm-A$3t 

+  lAC-S.C^  -  dz(Ci£j 

z  (B, BtCt-  EaC,- f  (BtBiCy&tt 

t  (  BiBjCi-BjCi  -B?Ct+  BiBjCtfSj 


fore, 

atlve. 


But  from  (Sx$)y“  ve  know  jfcfcat  tacBaO  for  any  values  of  "f  or  if.  There- 
,  except  for  the  case  Ik(bcC)s9,  the  vector  triple  product  is  not  aeaoci- 


The  triple  scalar  product  is  of  the  type: 


or  (Axt)-t 

Ve  can  immediately  racnpve  the  parentheses  since  (X*t)xJ  or  tx($**$)  is 
less.  Alsol*%6?  since  the  dot  product  is  commutative. 


1 


•>l»  wiiWMtw.  -  • 


*!W"  >  t« 


lot  ve  ean  staov  / .  flXC  ~  C'^XS  *8'CXA}  -to*: 


a*bxc, 


Ai  A,  Aj 

S,  8,  B, 


=  AX8*C*'£#C*^ 


1C  C%  C» .  ^ 

She  croea  product  is  distributive.  Thsfc  is,  jbd(S#)«  AX&MHC.  This  ease  can 
be  easily  demonstrated  by  expanding  the  two  sides  of  the  equation. 

matfksrm  nr  rotactso  cooRnnuiE  srazae 

The  following  diagram  illustrates  two  coordinate  systems,  a  g£  system 
and  an  ^  system,  in  which  the  system  has  been  rotated  with  respect  to  Xjj. 
system: 

fM 

Wi 


lh  each  of  these  systems  the  vectors  gj,  go,  -j  hare  unit  length* 

However,  if  we  represent  the  length  ofwevectfir  T£  u  the  ^  system,  then, 
as  is  indicated  in  the  diagram,  tbs  component  ofJl  parallel  to  Si has  a  mag¬ 
nitude  less  then  1,  and  similarly  for  the  component  of  parallel  to  ^* 
Thus,  if  we  are  dealing  with  two  coordinate  systems,  one  of  which  has  been 
rotated  with  reepect  to  the  other,  the  magnitudes  of  the  unit  vectors  in  one 
■yriem  will  not  ha  unity  if  described  in  the  other.  Let  ue  take  the  deriva¬ 
tive  of  a  vector  described  in  the  if.  system,  vector  X.  Where X  is: 

A«A,*(  *  A,*, 

9nJ 


sr: 


£ 


+  Aj ||  ^  +  (A,*?,  e*  Aj^) 


.  *  from  this  grouping  we  see  that  the  derivative  of  A  is  the  derivative  of  the 
ocapeaents  .i«n  eeoh  axle  still  in  the  direction  of  that  axis,  ths 
■  derivative  of  the  erne  themselves.  We  mow  with  to  obtain  an  eapreaaion  for 
determining  the  derivative  of  the  unit  vector*,  from  the  fallowing  diagram 


If  r]_  coincided  with  gi.  e.ndjr-}  coincided  with  gj,  then  after  a  rotation  through 
A$>  the  new  r^  and "r^.  say^'andTj',  would  be: 

t'-t+W  Zh<J  rhf=ti+Frf3 

If  A01s  very  small,  then  we  can  write  the  approximations : 

AT,  z -\r,\Mer>  and  Zr3-  JTJ/^a*; 

or: 

£ft  =  -AB%  a.nd  X?i=A^ 

Dividing  by/t  gives: 

and  $---??£ 


Which  in  the  limit  as  At-eo  is  exactly: 


This  gives  us  an  expression  for  the  derivatives  of  the  unit  vectors* 
However,  thfere  is  a  acre  convenient  method  for  determining  this.  If co  is 
the  angular  speed  of  rotation  then: 

AO  r  co At 

Define  the  vector  angular  velocity  as: 

60  =  CO  7^ 

that  is,  its  direction  i«  along  the  axLt  of  rotation.  The  cross  prodaot  of 
C5 and  the  various  vectors  is: 

.  ^  ’  r*  r, 

0  CO  o 
I  6  O 

r,  ft  r3  -  oSr, 

SCO 
O  0  i 

i3 


and  . 

tixr,  = 


But  By  definition: 


We  are  now  In  a  position  to  make  rules  concerning  the  derivative  of  • 
vector.  First  vs  have  the  derivative  of  the  scalar  components  of  the  vector— 
in  terms  of  the  rotating  coordinate  system.  To  tki£  "  must  add  the  cross 
product  of  the  angular  velocity  of  the  rotating  coordinau;  system  tines  the 
vector  itself,  as  defined  in  the  rotating  system.  This  gives  us  aa  aquation 
of  the  form: 


JL"t  X* 


oj  X  A 


where  the  superscripts  "g"  and  "r"  -n  >  it  the  derivatives  are  to  be  taken 
in  the  g;  and  •systems,  respectively; or =;t  the  superscripts  on  the  right,  j 
mean  that  .  the  ■  components  of  the  vector  are  to  be  expressed  in  the  indicatedo 
system. 


We  observed,  in  the  foregoing  analysis,  that  something  which  rotates  cam 
be  described  effectively  by  assigning  its  vector  direction  along  the  axis  of 
rotation.  If  the  rotation  i 3  positive,  it  would  advance  a  right-handed  screw. 
This  was  previously  defined  as  positive  in  the  coordinate  system.  Of  course  we 
can  immediately  define  a  negative  rotation  as  we  did  above. 


The  above  analysis  covers  all  of  the  essential  points  which  are  likely  to 
be  encountered  with  the  use  of  vector  analysis  is  the  study  of  missile  systems. 

The  algebra  and  the  study  of  vectors  covers  a  far  wider  field  than  we  have  done 
above.*  However,  vs  rarely  need  to  get  into  the  other  aspects.  It  is  likely 
that  ve  will  encounter  situations  where  we  need  a  greater  knowledge  of  vector 
analysis,  but  most  of  the  analysis  for  the  type  of  work  with  which  we  are  con¬ 
cerned  hare,  can  be  dona  by  observing  the  few  simple  rules  as  we  have  defined  them. 


EULER  ABQZJS  TRANSPOBMATIOH 


So  far  in  our  analysis  of  missile  systems  ve  have  only  described  an  equa¬ 
tion  for  the  forces.  In  this  equation  ve  needed  to  determine  three  types  ox 
farces.  The  gravitational  force  oould  be  immediately  written  set  •mg'fL.  Ve 
presume  that  ve  know  the  mass  of  the  missile;  certainly  g,  the  constant awoalar- 
atioa  due  to  gravity,  Is  kaceni,  so  the  gravitational  term  is  oceqpletely  defined. 
However,  we  do  not  know  the  thrust  as  such,  now.  Whether  cr  not  we  nhnoae  to 
put  the  thrust  In  a  thrust  system  with  finite  of  misalignment,  we  still 

auart  determine  the  thrust  in  the  body  system.  In  general,  ve  axe  given  some 
sort  of  curve  or  equation  which  will  describe  the  magnitude  of  this  thrust  vec¬ 
tor.  If  ve  put  the  thrust  in  the  body  system  (thereby  assuming  that  the 
between  the  thrust  coordinate  system  and  the  body  system  axe  aero)  ve  still 


•Fox  those  who  wish  to  pursue  the  subject  further,  there  are  any  eaoellent 
treatments,  for  example  see  Vector  and  Tensor  Analysis  by  Harry  Lass,  Nad ram 
HU1  Book  Oo.,  1??0.  ~ 
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sturt  achieve  some  method  of  ieoc:rf.bir!g  the  effect  of  the  thrust  on  the  body 
and  the  consequent  maser,  irara  that  the  alas  tie  win  undergo. 


She  missile  body  will  be  flying  an  arched  trajectory;  it  vill  be  mov- 
•  lag  frees  side  to  side  and  it  may  even,  be  rolling  back  and  forth,  or,  la  tom* 
cases,  spinning  around  Its  eva  axis*  Thus  ve  need  to  hare  some  sort  of  math¬ 
ematical  method  for  determining  the  orientation  of  tin  missile  body  at  any 
instant  of  time  with  respect  to  our  fixed  unchanging  axes  on  the  ground*  We 
mill  develop  the  attitude  angles  between  the  missile  and  the  ground  by  means 
of  what  is  known  as  an  Euler  angle  sequence.  Jfotbamaticajoy,  ve  assume  that 
ve  Smew  the  amount  of  angular  rotation  the  body  has  undergone  at  any  time* 
Then,  if  vs  are  able  to  completely  describe  all  subsequent  forces  which  causa 
the  missile  to  rotate  (thereby  saving  its  fixed  body  coordinates  with  it)  ve 
are  able  to  keep  track  of  the  change  in  these  angles,  and  ve  will  always  *nov 
the  orimtaticn  of  the  body  with  respect  to  the  ground* 

Standard  usage  throughout  this  country  in  the  various  airplane  and 
missile  industries  has  dictated  that  certain  symbols  are  to  be  used  for  cer¬ 
tain  angles*  Headers  who  have  used  spherical  coordinates  or  cylindrical 
coordinates  should  realize  at  this  point  that  the  coordinate  transformations 
ve  are  describing  are  net  the  same*  Although  these  coordinate  transformation* 
essentially  convert  a  rectangular  cooreurate  system  through  rotary  angles  to 
another  system  it  is  not  the  same  as  transforming  from  a  rectangular  system  to 
say,  a  spherical  one.  We  obtain  the  body  system  by  a  series  of  rotations;  a 
spherical  oocrdlrate  syjtsa  ia  obtained  by  projecting  various  vectors  onto 
planes,  the  components  sc  projected  determining  the  irigth,  direction,  mad 
size  of  the  saw  coordinate  syntax. 


Because  of  standard  usage  ve  will  make  the  following  definitions ;  if 
we  rotateabout  an  axf  ts  that  is  nuabersd  one,  aa  for  example  the  ~g±  axis,  or, 
say,  the  axle,  we  will  call  tilt  a  "roll"  angle.  Wheasrver  we  rotate  around 
a  number  two  axis,  cf  any  system,  vs  vill  call  this  a  "pitch"  angle.  Da  a  like 
manner,  any  rotation  about  an  axis  numbered  three,  ve  vill  call  a  "yav*  rat**- 
tico*  In  the  Euler  af.qr.^io-s  from  the  groxtd  system  tc  the  body  system  the  roll 
angle  is  defined  as  f,  ti>.  pitch  angle  io  defined  as  ©,  and  the  yav  angle  is 
defined  asT*  Geaaetriially,  vat  ve  era  trying  to  do  when  wa  transform  co¬ 
ordinates  is  the  following;  v«  start  at,  say,  the  ground  and  rotate  our  axis 
through  three  angles  so  that  they  are  all  pointing  in  the  same  direction*  as 
the  axes  of  the  ayatea  tc  which  ve  are  transforming.  (In  this^case  the  trass- 
formed  X  unit  vector  vill  point  in  the  seme  direction  as  the  unit  vector, 
the  trasSformei  *fg  unit  vector  in  the  safes  direction  as  the  b2  vector,  sod  the 
transformed  £  vector  will  be  parallsl  to  the  T-.  vector. )  Shis,  then,  is 
actually  wtat^ve  man  by  traasf cessation,  geometrically. 

Xa  physios,  as  generally  taught  in  school  systems,  trsaaf carnation  be¬ 
tween  QKXCdinate  system*  Is  performed  by  a  roll,  pitch,  roll  sequence.  Cm 
reason  for  this  particular  type  of  transformation,  X  do  sot  knew.  3fc  our 
atek  ve  will  sot  restrict  ourselves  to  any  particular  type  cf  transfers*  tl  on. 
Qaa  sequence  of  angles  is  as  good  as  another  sequence  with  one  restriction* 

For  example,  suppose  we  decide  to  go  through  a  roll,  pitch,  yav  sequence*  Than 
the  restriction  cm.  this  soqvtence'  is,  that  the  second  eagle  (in  this  oasa,  ©), 
can  never  be  Mar  90*  The  reason  fax  this  is  that  the  second  angle  appears 
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In  the  final  transformation  as  a  tangent.  If  the  second  angle  nears  90°,  than 
the  tangent  of  this  angle  approaches  infinity.  Other  than  the  restriction  on  „ 
the  second  angle  we  can  take  such  sequences  as  0,  <p  which  is  a  pitch,  yaw, 
roll  sequence.  We  can  take  a  sequence  such  as  0,  0,  *)P*.  This  is  a  pitch,  roll, 
yaw  sequence,  etc.  No  matter  which  type  of  sequence  we  care  to  use,  we  must 
adhere  to  the  definitions  of  pitch,  yaw  and  roll,  that  is,  the  pitch  is  around 
a  number  two  axis,  a  roll  is  around  a  number  one  axis,  and  yaw  is  around  a  num¬ 
ber  three  axis.  The  methods  used  for  all  transformations  between  the  ground 
and  body  by  an  Euler  sequence,  are  the  same.  Certainly,  the  method  of  solring 
the  equations  will  be  similar. 

For  the  purposes  of  demonstration  of  the  methods  employed  to  transform 
the  ground  to  body  we  will  use  a  ©,  ‘Y-,  0,  sequence;_J&at  is,  a  pitchy  yaw,  roll 
sequence.  Th ' s  means  that  ve  first  must  rotate  the  system  around  gg.  From 
this  we  will  define  the  rectors  of  the  subsequent  coordinate  system  in  terms  of 
the  gp  vectors.  Ve  will  write  down  the  vectors  in  the  following  manner:  "gpi  means 
the  vector  after  the  first  rotation,  "g^  will  be  the  "gg  vector  after  the  first 
rotation,  and  will  be  the"gh  vector  after  the  first  rotation.  Similarly 
will  he  the  %  vector  after  the  second  rotation,  etc.  Also  during  these  trana- 
fonnationSjVe  must  keep  track  of  the  angular  velocities  so  that  we  can  determine 
the  angular  velocities  of  the  missile  in  terms  of  the  missile  coordinate  system- 
chat  is,  the  bedy  system.  The  following  diagram  demonstrates  how  we  obtain  the 
coordinate  system  after  tho  first  rotation,  vhich  in  this  case,  is  a  pitch  ro¬ 
tation. 


Wa  can  write  the  g^s  in  terms  of  the  (£^3  Uf 


fu  *3,C*S6-gsJin* 

5.*?* 
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Having  obtained  the  coordinate  system  which  vs  trill  call  the  system  after 
the  first  rotation  we  will  now  rotate  through  the  second  Siler  angle,  'f',  a 
yaw  rotation, '  in  the  "g^^syste^.  We  must  rotate  about'  a  number  three  axis, 
so  we  will  rotate  about  831*  The  following  diagram  Illustrates  the  derivation 
of  the  coordinate  system  after  the  second  rotation: 


tzr  C6S  ? 

■A  ' 

9j  l  z 

wc.  m.*Y23i 


In  order  to  transform  the  ground  system  axes  so  that  the  corresponding 
numbered  unit  vectors  are  parallel  to  the  corresponding  numbered  unit  vectors 
of  the  body  axes,  it  is  necessary  to  use  three  distinct  angular  rotations.  We 
'  have  used  two  angular  rotations  so  that  one  more  rotation  is  required  to  align 
.  our’  ground  system  with  the  body  system.  The  final  rotation  is  a  0  or  roll  angle. 
This  means  we  must  rotate  about  *g-2  that  is,  the  first  ^  vector  after  the  second 
rotation.  Since  we  wil.'  ov  be  aligned  in  the  ’^"System:  ou5T third  rotation  is: 


J?i-  3ig 

-  Qzz  cos<f>  +%iSin<l> 

%z  -%zM<P+’3mc*s$ 


AUG-  VEL.s  <pb, 


We  now  have  expressions  relating,  by  a  series  of  rotations,  the  ground  system 
to  the  body  system.  Also  we  have  kept  track  of  the  angular  velocities  for  each 
rotation  which  we  performed. 
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The  form  of  the  Euler  angle  relationship*  i*  awkward  to  use,  so  a*  a  cos- 
pact  tool  we  will  introduce  the  transformation  matrix  relating  the  ground  cod 
body  vector*.  In  matrix  station,  the  "$£  vector*  are  related  to  the'gj^  vec¬ 
tors  's&t 


In  a  like  aaaner  vs  can  represent  the  vectors  of  the  g^'a  after  the  second 
rotation  as; 


COSf"  smf 

-smf 
c  o 


and  the vectors'ia  terras  of  the  gj  vectors  after  the  second  wst®'  m  urn 


^  M 

h 

0 

be 

i 

f  c- 

cosj> 

W 

1° 

-sir>$ 

I  Y&2/ 


If  we  represent  these  matvix  equation*  in  a  nor*  coolant  notation  w  can 
write  the  following  equations: 


(Jt}=  M  (©)(&) 
Cic^’M  (r$<) 
(b i)--M  (<PMu) 


Then  substituting  in 


the  equation  for 


sad 


we  obtain  the  following  *$**» 
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iion: 


(bi)=  tf 


Rathar  than  indicate  the  multiplication,  of  these  three  matrices  together  ve 
will  give  the  final  result.  That  is: 


f~&\  fcosy&se  sin  y 

uJ 


■  CilS'f’SlllB 


\  h\ 


\ 


bj 


7 


\ 


cos <HinrceSMinisth& ~  cosfcesf  css 4 HflfSitl0+Sin4 
SHtiSinfeue+ttSfar,*  -sinfasf  -sir\$-*Mi'Sin9+£c£lfeeS^' 


3z 

V  */ 


Derivation  of  the  Suler  Angles 

■  The  above  matrix  identity  gives  us  immediately,  presuming  we  know  the  site 
of  the  the  three  angles,  the  relationship  between  the  ground  and  the  body  coor¬ 
dinate  systems.  As  everyone  is  probably  already  aware,  every  element  of  the 
matrix  on  the  left  side  of  the  equ-.-I  sign  must  equal  the  element  corresponding  in 
position  in  the  matrix  on  the  right  3ide.  For  example  "b^ ,  from'  ;  the  product 
‘of  the  two  matrices  on  the  right  side,  iabuq: 

b,  =  COSYCQS&g,  Jr$ir>yfi - COS^SthB-^ 


Similarly,  we  can  find  expressions  for  bg  and  b^.  Thus  if  we  can  write  a 
vector  in  terms  of  the  body  system,  we  can  then  employ  this  transformation  ma¬ 
trix  and  express  it  in  terms  of  our  fixed  ground  frame.  However,  we  have  in¬ 
troduced  three  more  unknown  quantities;  'fi  Q,  0.  Therefore,  further  analysis 
will  require  that  we  find  methods  and  relationships  in  order  to  solve  for  these 
angles* 

The.  following  angular  velocity  notation  for  the  bjdy  system  is  a  standard 
adopted  by  most  of  industry.  The  angular  rate  around  b]_,  the  roll  rate,  is 
designated  by  P,  the  angular  rate  about  Dg,  the  pitch  rate,  is  represented  by 
Q,  and  the  yaw  rate  by  R.  The  angular  velocity  vector  can  be  considered  as  a 
vector  in  space.  Ve  have  two  coordinate  systems^  vail -defined,  with  the  trans¬ 
formation  between  them.  We  can  write  the  angular  velocity  in  terms  of  either 
coordinate  system,  but,  because  of  ease  of  computation, which  ve  will  see  later, 
we  will  solve  for  the  angular  velocity  with  components  in  the  body  system. 

Since  we  kept  track  of  the  angular  velocities  during  our  transformation,  we  can 
write: 

*  Note  •  that  the  transformation  matrices  as  defined  can  be  utilised  for  other  Euler 
sequences.  For  example,  a  sequence  would  be  written  as 

M(0)  dj)*  angular  velocities  would  be  written,  .however  as: ^3,  06u, 

and  dbg. 
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pt|+  Q  V  -  *0%  +  *%+Y«xm  £  V  t*i*  9%*Teoa  9% 


V*  in  new  able  to  solve  for  the  angular  rates,  as  rogroaeated  la  the  body  sys¬ 
tem,  in  terns  of  trigonometric  function*  and  rate*  of  the  Solar  agios,  th§% 
selves.  taking  the  dot  or  scalar  product  on  each  aide  of  the  equation  within 
give*: 


P  »  0  +  0  ggS  +  y®1»  ©  g| .  *b,+  'fcos  0 


In  addition  ve  can  find  Q  by  taking  the  scalar  product  on  sesh  aids  of  the 
equal  sign  with  "Sg ,  and  ve  can  find  R  by  taking  the  scalar  product  with  on 
each  aide  of  the  equal  sign.  This  gives  ua  the  following  aeries  of  three 
equations: 

P  «  0  +  ©"gj.  %t  *  ^a±a  9  "gf  ."bj  ♦  ‘f’teos  O'g^'hj 

Q  »  0  g^.  *^+  ywln  ^eoa  0^.*^ 

R  *  0  Y*ln  Ogj .  "tyt  'fee*  •VI 


Proa  the  transfonaati'w  matrix  which  ve  derived  ve  can  iamadietely  find  the 
dot  products  of  the  vectors  necessary  to  solve  these  equations.  Substituting 
in  the  appropriate  trigonometric  functions  for  all  of  the  dot  products  gives 
us  the  follovii^  expressions  for  P,  Q  end  R: 

P-  4> 

Qsh  CAS+COS  Y~  +  Y  Strf&JHlj 

<3-  ec6s$c0$Y+ym$ 

Rs  -&j//rfa>JY+  Y&jJO&fmYfmj  + $co$4>sirf* 

,  +Yc&bcAsf~Vcas*f&j*ttfYfin0‘ 

Early  in  these  series  ve  mentioned  that  ve  would  require  six  degrees  of 
freedom.  So  far  ve  have  investigated  throe  of  them.  These  art  the  throe  com¬ 
ponents  of  the  force  vector.  The  other  three  degrees  of  freedom  arise  from 
the  three  components  of  the  momentum  equations.  The  Momentum  equations  give 
xis  an  expression  for  obtaining  P,  ft,  and  R.  Therefore  if  ve  assume  ve  can 
find  P,  ft,  and  R  ve  can  use  the  above  equations  for  finding  0,  f,  and  0* 
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jQat  matrix  notation  the  &tore  can  be  written  eat 


/P\  !\  sm V  o  \l  j> 

|  Q  Jr  jo  Ctefcesf  Sin4  I  & 

\rJ  \0  V 

lh  compact  notation,  represent  the  above  equation  by: 

/!  -  M  R 

At  this  point  we  should  digress  scsnwbat  to  explain  the  tam  Orthogonal " 
matrix.  In  order  to  solve  the  shove  matrix  equation  for  °0,  6,  and  "Y  we  must 
find  the  Inverse  of  tint  matrix  which  relates  P,  Q  and  R  to  the  derivatives. 
Tte  Inverse  is  found  from  the  equation 


R  -  vT^fl 

*  *  t 

which  gives  us  an  expression  for  'Y,  Q,  and  0. 

In  the  analysis  which  we  did  for  the  transformation  from  ground  to  body 
coordinate  systems,  we  obtained  the  transformation  matrix  whloh  related  the 
body  and  the  ground  systems.  This  transformation  matrix  possesses  a  unique 
property*  because  of  Its  nature,  the  Inverse  of  this  matrix  Is  equal  to  its 
transpose.  If  we  write  the  relationship  between  the  ground  and  body  systems 
as: 

then 

GsK'MqQB. 

However,  for  this  matrix  we  can  prove  that 

When  we'  can  do  tbit*  the  matrlx/  ln  thls  case:  the  transformation  matrix,  is 
call  ad  orthogonal. 

T ho  matrix  relating  the  angular  rates  in  tbs  body  system,  P,  Q,  R,  with 
the  derivatives  of  the  Buler  angles  l.v  not  an  orthogonal  matrix.  Therefore, 
we  asm  cospelled  to  derive  its  inverse.  Without  goli«  Into  the  actual  com¬ 
putations  of  this,  the  Inverse  matrix  1st 
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-smycos* 

cos$ 

MfGUy 


Sinystn#' 

~sm$ 

cuY 


Ws  in  now  In  a  position  to  know  at  any  Instant  of  time  ths  anglss  yst  Qt 
$  ,  merely  ty  Integrating  their  derivatives  w'.th  reopeot  to  tin*.  However,  at 
this  point  we  have  merely  replaced  the  three  unknowns  *y»  0,  0  ty  three  other 
unknowns  P,  Q  and  R. 

In  tody  coordinates  we  can  write  the  force  vector  as: 

We  can  solve  tor  each  of  the  components  of  the  force  in  the  body  systssi  ty  the 
simple  expedient  of  taklsg  the  dot  product  with  "£j>,  andt*,  successively. 
Doing  this  we  obtain  the  follcwlig  thiec  aquations: 

F,  -  -  ■ %•  b,  +  %  •£,  + 

Fix  -rw&'biL  i'T^bx  +%•  %. 

fS-  -**3$*%  +%'%+  F r*% 

Since  we  already  hsnre  the  transfoxvatlon  matrix  between  the  ground  and  body 
system  we  can  write  down  the  trigonometric  functions  which  relate  the  gravi¬ 
tational  force,  f^p,  to  the  body  system.  This  gives  us  the  series  of  equations: 

F{  *  rngcos‘tJM*  +  fnt6)  +Fr'*t 

ft  m  -frtjfcos+jin+*m0tJin+C6S*)  t  +  %*%* 
f9z  -ntp  (CMfast-iir^sin  % •% 


tBX  HUD  SX9TKK 

As  yet  we  do  not  know  what  the  transformation  is  between  the  —  rmljimmln 
foroe  vector  and  the  body  system,  tfs  represented  it  above  ty  the  sealer  product 


S3 


i 
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Tjetween  each  unit  vector  In  the  body  system,  and  the  aerodynamic  force  rector, 
F».  As  stated  before,  the  ^aerodynamics  are  represented  In  the  vlnd  system, 
vnich  la  defined  so  that  "v.  always  points  parallel  to  the  relative  wind.  We 
will  define  the  total  relative  wind  vector  as  Wy.  The  vector  1«  represented  as: 

where  the  "Ww"  Which  is  not  a  vector  Is  the  magnitude  of  the  relative  wind. 

The  aerodynamic  forces  arise  entirely  from  the  effect  of  the  relative  wind 
age.1  net  the  miseile  body.  Therefore  we  need  to  represent  the  relative  wind  in 
a  system  which  we  can  define.  The  force  acting  an  the  missile  will  bo  a  force 
parallel  to  the  dl&octlon  of  the  relative  wind.  In  order  to  perform  this  ana¬ 
lysis  we  will  transform  the  body  system  Into  the  wind  system  by  rotating  through 
two  angles,  celled  the  angles  of  attack.  These  angles  are  o c  and  .  o<  is 

called  the  pitch  angle  of  attack  or,  often,  the  angle  of  attack.  la  called 
the  yaw  angle  of  attack  or  the  sideslip  angle.  Collectively,  we  call  oc  and 
the  aisles  of  attack. 

Previously  when  we  were  transforming  from  the  ground  to  the  body  system, 

It  was  stated  that  three  Euler  angles  were  required  In  order  to  align  the  three 
ground  axle  vectors  parallel  to  the  corresponding  body  vectors.  In  the  trans¬ 
formation  from  the  body  system  to  the  wind  system  we  are  Interested  in  aligning 
the  ox  rector  with  the  ^  vector  only.  To  do  this,  we  only  need  to  rotate 
through  two  angles.  The  body  two  and  body  three  vectors,  of  course,  will  be  ro¬ 
tated  at  the  same  time,  but,  because  we  are  Interested  in  aligning  only  one 
vector,  wj, will  call  thatfc,  and*^  u^lt  vectors,  those  vectors  which  result 
from  the  bg  and  03  vectors  after  rotating  the  body  system  through  the  two  angles 
of  attack.  This  simplifies  the  derivation  a  great  deal,  since  we  needn't  define 
a  third  angle. 

The  method  of  transformation  Is  similar  to  the  on s  used  before.  First  we 
will  rotate  through  s  pitch  angle  of  attack,  that  is,  PC  ,  and  then  through  the 
year  angle  of  attack,  f  .  Again  we  will  observe  cur  criterion  for  positive  ro¬ 
tations.  The  following  diagrams  detenslns  the  trenafoxmations  required  to  go 
from  the  body  system  to  the  wind  syetsm: 
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"b,,  -~b3stnoc 

hsr  'StSMoc  h£$co6ac 
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%\CO$/5+b2iSifl/ 

bsi 


Tram  these  expressions  we  obtain  a  transformation  Matrix  between  the  body  and 
the  wind  systenu  Ibis  Matrix  (which  la  orthogonal)  la: 


I  C&SOL  CAS  0  Sinf  -Stnou; af/*\  /  % 

-cosxsmfl  CASfi  Jtnocsw 
\  sincx  o  C6soL 


Utilising  the  transformation  Matrix  which  we  found  during  this  rotation, 
we  can  now  represent  the  aarodynanlc  forces,  which  are  defined  in  the  wind 
system,  in  teres  of  the  body  coordinate  vectors.  Thus  we  can  write  the  fallow¬ 
ing  equation: 


fA  ~  Puj,  +  Su/g-b  L&r-CbitYhz+Nbj 


The  teres  D,  8  end  L  in  the  wind  eyetea  are  given  the  following  nanee:  0  is 
the  drag,  3  la  the  side  force,  and  L  Is  the  lift.  These  are  the  names 
which  ws  apply  to  tbs  aerodynes lo  forces  in  the  specific  direction  Shown  stove, 
when  they  axe  defined  In  the  wind  eye tea.  The  teres  In  the  body  syeten  are: 

C,  the  chord  foroe,  X,  the  yew  force,  end  I,  the  noxnal  fores.  These  are  the 
naaee  given  to  the  speclflo  eerodynenlc  forces  defined,  jes  stove,  lathe  body 
coordinate  eyetcau  If  we  now  take  the  dot  product  vith*e»  ,T>p,  and  03,  tr  ' 
eucceeslvely,  we  obtain  expressions  for  the  chord  foroe,  the  yaw  fores,  and 
the  normal  foxes  in  terse  of  the  drag,  side  fores,  and  lift.  The  dot  protests, 
la  terns  of  tbs  trigone— trie  functions  of  the  angles  of  attack:,  sea  bn  innedlately 
found  from  the  trvnsfoneeti  an  netrlx  Which  wo  >ist  derived. 
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Derivation  of  Angles  of  Attack 


So  far  w  still  sees  to  be  behind.  We  sey  assume  that  the  drag,  side 
force,  and  lift  ere  functions  cr  curves  which  ve  know,  hut  ve  have  introduced 
two  unknown  quantities,  the  angles  of  attack.  However,  ve  can  detenri.ro  these 
angles  of  attack  rather  ei*ply  by  means  of  the  following  derivation: 

We  can  represent  the  relative  wind  vector  aa 
^  =r  V^y^i  +  + 

whexe  Vn ,  Wg,  and  V*  are  the  components  of  the  relative  wind  defined  in  the 
body  aystae.  Taking  the  scalar  product  of  the  relative  wind  with  ^  gives: 

K- Hb,&,  +  ^ h'v;  +  tiJb/Z, 


which,  when  ve  utilise  the  transformation  matrix  Just  defined,  gives: 


Vfc 

W, 

From  thia  wa  define  aln  (-  0(  )  and  co« ( - 0(  )>  whit*  are: 

sm(-*)z-smoc  s 

CO$(-<x)zCjQJ<X-  ^ 

It  va  now  Multiply  the  firac  of  the  prerioua  equatiana  by  tin  jff,  the  aaocad 
locr  ooa  ,  *?&  add,  we  obtain 

ty  co$0(Sirtfccjfi+ Vk  sirfy-  tysmxsmficaifl 

_ 0  -  -K  casQLS\nPo#fi+ W? 

Wfe 


If  w*  multiply  tha  firat  aquation  tgr  eoa  fl  ,  tha  aaoend  hy  ala  ,  and  «V 
traet,  w  obtain  ^ 

W*0(?y-  casocoxf*^  t nk* ~ 

WwCa3/f?  ^casK-tfjJW  t  , 

*w^  «  <$)  - 

•'•  WV"  $r 
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Thus,  vr  have  obtained  exact  expreeaiona  for  the  angles  of  attack,  C*  ,  ^8,  1" 
terms  of  the  components  of  the  relative  wind  in  the  body  system.  The  angles 
of  attack  are  no  longer  unknowns;  the  terms  of  the  relative  vlnd  ere  the  new  un¬ 
knowns. 
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The  only  force  which  has  not  been  defined  In  the  body  system  at  th'j  tins® 
is  the  thrust.  Previously  we  said  ve  were  defining  the  thrust  so  that  it  was 
parallel  with  ap,  vector.  It  is  necessary  for  the  person  who  is  doing  a  study 
with  thA  missile  to  decide  his  own  thrust  misalignment,  ibis  means  that  the 
thrust  vector  (which  always  acts  in  die  plant  of  the  -exit  noasle,  not  at  the 
throat)  cm  hare  a  directional  error  (by  error  we  mean  it  does  not  act  along 
the  missile  center  line)  as  well  as  an  assumption  that  the  thrust  vector  does 
not  act  at  the  center  of  the  exit  nosale.  Thus  ve  can  create  an  error  in  posi¬ 
tion  as  In  the  following  diagram. 


as  well  aa  a  directional 


tour  AmicAri»* 
«r  thh  -nUtonT  Were* 


error  as  shown  in  the  following  diagram: 


The  analyst's  discretion  determines  how  large  the  thrust  misalignment 
angles  win  be.*  in  general,  thrust  misalignment  angles  an*,  at  most,  a  few 
degrees.  Ouce  ve  choose  the  nlsallgnaent  angles,  ws  can  create  a  transforma¬ 
tion  matrix,  aa  we  did  before,  between  the  thrust  axis  ayst**  and  the  body 
system  —by  transforming  the  bj  vector  into  the  vector.  A*  vlth  the  aero- 
dynamics  utilising  the  angles  of  attack,  we  need  only  define  two  ^ioalignaent 
angles.  If  ve  cell  T p  the  pitch  and  y  the  yew  misalignment  nngles,  and 
assume  a  transformation  Tree  the  body  to  the  thrust  system  by  means  of  a  f  p. 
Tv  sequence,  the  transformation  aatri.c  is  precisely  like  that  for  the  wind 
ayvUm  with  ot  replaced  by  f  p  end  replaced  by  J"y.  'rhui! 
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•  In  some  cases  It  is  possible  that  sons  sort  of  thrust  misalignment  angles 
will  be  furnished,  but  in  general  we  will  investigate  a  spectrum  of  values 
between  some  reasonable  limits  chosen  by  the  analyst. 
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If  the  thrust  mi ^alignment  angle*  were  zero,  of  coirse,  this  implies  (the 
more  usu*1  case)  that  the  thrust  ects  along  the  missile  center  line  and  is 
parallel  to  the  vector.  In  any  case,  because  of  the  way  missiles  are 
manufactured,  we  can  assume  tv*t  the  thrust  mlsallganent  angles  are  very 
small,  and  utilize  the  "small  angle"  approximation.  This  gives  the  transfor¬ 
mation  matrix  the  form: 


fy 

I 

o 


-r) 

0 


where  we  aesumed 


&fy*0. 


If  we  presiaMi  that  we  h#  m  same  sort  of  curve  or  function  wbloh  repre¬ 
sents  the  thrust  characteristics  of  the  actual  missile  with  which  we  way  he 
dealing,  we  can  represent  the  thrust  as: 


FT-rfi-  Tibi+Tibt+Tjbj 


Taking  the  scalar  product  with  th, sequentially,  wo  obtain  three  equa¬ 
tions  for  the  magnitudes  of  the  thru-i  parallel  to  the  body  coordinate  axis: 

77  =  Tf,‘%  -  T 
Te  =  Tf,’%=Trr 

Ts--Tf,'%  =  -Tr, 


t  raaoncs  a kd  rsutivb  void 

At  this  point  many  unknowns  have  been  solved,  but  it  appears  that  as  fast 
as  we  solve  thee  we  find  new  ones  to  taka  their  places.  Oerta-nly  this  Is 
true  with  the  relative  vind.  Before  we  define  the  relative  wind,  which  la  a 
combination  of  the  missile  velocity  with  rsepect  to  the  inertial  flees  and  say 
external  winds,  ws  will  investigate  the  missile  velocity  with  reepecrt  to  the 
ground  system. 
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f.i*K.'  TranolatJ  on  /LhvJ  e  iv-tlonu 

The  differential  equations  describing  the  motion  of  a  missile  system  must 
always  he  referred  to  the  ground  or  inertial  reference  frame  for  solution.  Bile 
is  not  meant  to  imply  that  every  computation  is  done  In  the  ground  system  ■—  to 
the  contrary,  most  computations  are  performed  elsewhere.  It  does  mean  that  the 
basic  equations  of  motion  must  he  referred  to  the  inertial  system  and  computa¬ 
tions  in  other  systems  related  to  the  inertial  frame  by  appropriate  trensfomna- 
tion. 

To  develop  the  transformations  between  the  inertial  frame  and  rotating 
framt-  we  ere  at  liberty  to  choose  a  vector,  assume  it  is  in  space,  and  define 
it  ii_  terns  of  its  components  in  any  cooruinate  system  which  we  desire.  But 
once  we  have  chosen  such  a  vector,  subsequent  differentiation  and  integration 
must  involve  the  motions  of  idle  rotating  frame. 

When  we  describe  a  missile  system,  ve  are  never  likely  to  need  range  ex¬ 
pressed  in  body  coordinates,  but  we  will  need  the  range  in  the  ground  frame. 

We  do  need  the  velocity  in  body  coordinates  because  the  aerodynamics  are  based 
on  missile  velocity  through  air,  so,  in  general,  for  our  kinematic#  we  wish  to 
define  motion  by  means  of  the  velocity  vector.  Thus: 


dr 

VA+ VV  f 14& =u£,t  v$e  +  iJb3 


If  we  know  the  velocities  in  the  body  system  then,  Simply  finding  the  dot 
products  with  "gp  "gg,  arid  re«  actively,  we  can  obtain  the  velocities  in 
the  ground  system.  The  terms  u,  v,  w  are  standard  accepted  symbols  for  the 
velocities  in  terns  of  the  indicated  body  coordinate  axcc.  The  notation  above, 
that  is  the  superscripts  ”g"  and  "b",  indicate  in  which  coordinate  system  these 
vectors  are  defined.  (Note  that  it  is  not  always  necessary  to  state  in  which 
coordinate  system  a  vector  is  defined). 

According  to  the  above  equations  the  velocity  vector  In  the  ground  system 
is  directly  equal  to  the  velocity  .actor  in  the  body  system.  (We  could  have 
chosen  the  range  and  represented  it  in  the  same  way,  or  we  could  hwe  chosen 
the  acceleration).  In  order  to  derive  the  kinematics,  we  will  differentiate 
the  velocity  vector  in  the  body  system  for  the  acceleration,  integrate  it  in 
the  body  system,  and  then  transform  It  to  the  ground  system. 

We  previously  found  an  expression  relating  the  forces  to  the  accelerations. 
Above  we  achieved  an  expression  relating  the  velocities  in  the  missile  to  the 
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velocities  in  the  ground.  We  will  now  derive  the  acceleration®  in  the  body 
system.  Newton's  lavs  of  motion  are  valid  only  in  an  inertial  system,  thus  ve 
nruet  differentiate  the  expression  for  the  velocity  in  the  Inertial  system* 

That  is: 

s<Ly9,  14-\/i 

a£tv 


It  is  'immediately  apparent  that  the  expression  on  the  right  side  of  the  equal 
sign  Is  mixed  and  we  can't  work  with  it  as  it  is.  We  define  the  derivative 
with  respect  to  time  in  the  inertial  or  ground  system  as  the  operator: 


j 

i 

so  that 


It  is  now  apparent  *hy  we  vent  to  such  trouble,  previously,  to  illustrate  the 
equivalence  between  differentiation  and  cross  product  multiplication  in  a  rotating 
coordinate  system. 

Completing  the  above  differentiation  gives: 

‘4«A  +•  vri)  z-r  krbj)  -Mb,  b2  b3 

P  Q  R 
a  sr  w 

or 

f=  tj]?'  +  \r%t  +■  i-(Qur-R V^S)  t  (Hu  ■PurfSg 

+(fY-Qu$b3 


3<L 

2t 


bj  — 

&£.  +  co  y 
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anti  finally: 


“tfs  (u-PQus-PV^b,  -h  (^+Ru-P<Ar)b2  +  (&  +  Pf-Q&yhj 


* 

Solving  the  force  equation: 

— mJ *A  5  *”*7^ 

=  77)  V 
or 

JxZ.  +  El%  ±  &  jZi~(k+QM'R\r)b,+(fr+Ru-Pu/)k>& 

W  1  Jt}  *  7*1  J 

+(uf+P\T-Qu.)bj 


This  supplies  ua  with  the  three  Jtranslatienal  equations  of  notion  when 
vs  form  the  scalar  produets  vith*t>i,  b2,  and'%q  respectively.  Thus: 

d=  &r-Qur+RY 

y*i 

sr=Jg--Ru.+Pur 

u/s  §--P\ri-QU' 


J2 


The  Rotational  Accelerations 

In  the  last  /section  we  developed  three  very  pretty  equations  for  the  trans¬ 
national  accelerations  in  the  body  system  However  ve  still  have  the  three  un¬ 
knowns.  P,  Q,  and  B  which  must  be  determined.  7,  Q,  and  B  are  also  used  for 
computing  'f,  Q,  and  0.  These  ter®^  which  arise  from  the  rotation  of  the  body 
system,  require  three  additional  independent  equations  for  their  solution.  We 
can  find  three  such  equation?  because  the  rotational  accelerations  are  inde¬ 
pendent  of  the  translational  accelerations.  As  was  done,  in  the  last  section, 
we  find  the  rotational  accelerations  by  solving  Newton's  second  law  of  notion, 
in  terms,  of  the  rotational  accelerations. 


We  can  write  the  equation  for  the  torque  as 


f=FSJ 


where  T’ is  the  torque^d/dt  means  the  derivative  Is  taken  in  the  ground  coor¬ 
dinate  system,  and  (&)  is  che  angular  momentum,  where  I  is  the  ament  of  inertia 
tensor,  and  to  is  the  angular  velocity. 


This  rotational  equation  exactly  parallels  the  translational  equation  pre¬ 
viously  defined.*  X,  the  moment  of  inertia,  corresponds  to  the  mass  of  the 
translational  equation,  and  to ,  the  rotational  velocity,  corresponds  to  the 
translational  velocity.  As  pointed  out  in  our  previous  dlscusslo/cu  in  order  to 
take  the  derivative  in  the  ground  system  ve  must  replace  d/dt  of  the  ground  system 
by  its  equivalent  operator  in  the  body  system.  Thus: 


which  gives: 


’Z(i3. 


Note  that  in  the  abovT  equation  ve  have  the  teraG$x(iai).  This  cannot  be  written 
aa  I (Sbx&j)  (which  is  equal  to  zero)  even  though  the  inertia  might  be  composed  of 
constant  termsj  because  X  is  a  tensor.  In  general,  "X"  Is  represented  as  a  matrix 
(or  tensor)  containing  nine  terms.  Thus: 


/I„  -  hi 
”Tu  Til 
rlji  -In 


•Its 

I$3 


*  The  translational  equation  is:  ifcS/at  W),  that  is,  the  tine  rate  of  change 
of  the  linear  momentum  is  equal  to  the  sum  of  the  external  forces. 
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Bowever,  in  erary  rotating  physical  system/  It  cap  be  prwod  that  there  exist 
three  specific  sixes  in  which  all  terns  except  those  cat  the  nain  diagonal  equal 
serb.  Thm  m  can  obtain  the  following  matrix  using  a  eultahle  tranifa^tidh•: 


0. 

fid 


0 


2a  wldeh, if  rotation  occurred  about  these  axes,  the  ''cross  product  of  Inertia" 
fraud  be  sero.  Theee  axes  are  called  the  "principal  axes  of  inertia*.  Xt  is 
always  possible  to  find  such  axes. .  Generally  with  a  missile  system,  these  axes 
■not  be  found  by  experiment.  Sheerer,  if  we  take  a  hasogaasous  symmetrical 
body,  than  we  oan  locate  these  principal  axes  of  inertia  by  Inspection.  If  we 
assume  that  our  else  lie  system  is  a  eyamatricsl,  honogeneous  body,  then  the 
principal  axes  of  inertia  would  lie  along  the  sans  axse  ag  the  body  system 
axes  —  one  of  the  primary  reasons  vs  chose  to  work  in  this  system.  Of  course, 
the  missile  is  not  a  homogeneous  body#  but  for  stability  reasons.  It  Is  generally 
rery  dose  to  being  one.  Almost  always#  therefore#  we  can  assume,  quite  aecur 
retell  that  the  principal  axse  of  inertia  and  the  body  axes  coincide.  If  it 
weft  neeessoftr  to  be  more  exact  and  actually  use  the  oap&ste  font  of  the  moment 
’  of  inertia  tensor,  the  easiest  method  for.  analysis  would  be  to  locate  the  prin- 
*  clpal  axes  of  inertia,  define  the  unit  rector*  parallel  to  them,  and  solve  our 
t  rotational  equations  in  that  system.  This  would  necessitate  a  transformation 
from  the  system  in  which  the  principal  axes  of  inertia  are  located,  to  the  vody 
system.  This  involves  a  great  deal  of  extra  work  (in  the  ernes  way  that  trans¬ 
forming  from  a  nonrotating  system  at  the  center  of  the  earth  to  a  rotating  one 
on  the  surface  of  the  earth  requires  extra  work)  and  so  we  will  Just  mention 
this  fact  but  not  go  Into  details  of  calculations.  If  it  ever  becomes  necessary 
for  the  analyst  to  deal  with  the  general  moment  of  Inertia  tensor,  he  must  use 
the  seme  concepts  that  we  have  been  developing  throughout  this  paper. 


Assuming,  ngjr,  that  the  principal  axes  of  inertia  and  tha  body  axes  coin¬ 
cide,  tha  tern  Ico  is* 

/l«  0  °\lp}\  I3* Pr \ 

X  o)  r  ©  Xu  0  I  Qi*  Is  I  liiOtH  I 
\o  o  XmJ  \  Rht  I  \ln  fwsj 

furthermore,  if  the  missile  is  shaped,  roughly,  like  a  cigar,  rotation  about 
i#  is  Indistinguishable  from  rotation  about “bj.  Shat  1st 
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the  moment  of  Inertia  about  b-j,  and  that  about  *>2,  are  sjjbstactially  the  same. 
Therefore,  we  can  eet  Ig2  »  I33  and  this  ajives  us  for  XcO: 

•  i^Pb^  +  IpQbg  +  ipa3 

where  we  have  put  | 

Xgg  -  I33  -  Xp. 

We  .can  now  write  the  equation  for  the  torque  as: 


Vztf 

q  n ■!  • 

(i„p  w 

treating  the  moments  of  inextia  as  we  did  the  mass  (that  is,  we  will  extract  the 
Instantaneous  value  of  the  variable  when  we  use  it)  we  obtain: 

T  -  Ill*!  ♦  [jW+HP(Xir*P)J  *2  ♦  [W?<i(3>-Iii)]^3 

If  we  define  -  Iji  •  le,  then  the  torque  la: 

T -  Ill^i  ♦  +  &U)%  *  (Ip»  -  *42,5*3 

Another  name  for  the  torque  is  "moment  of  momentum".  A  physical  law  c 
states  "the  time  rate  of  change  of  the  angular  momentum  (the  moment  of  momentum) 
is  equal  to  the  sum  of  the  external  momenta”.  In  equation  fora  this  Is: 


T  -**/«  (iix^i+Ip^+IpBbj)  + 


b, 

P 


Jd/dt  (ICO)  -ZN 


or  component  wise 

T  "  #*/*t  (l<0)  -  Z*M  -  ♦  2Kb?  ♦  ZHb3 

where  21, 2M,  ZN  are  staodard  notation  for  the  external  moment  compocUnts  in  the 
Indicated  body  directions.  Taking  the  dot  product  of  the  components  of  the  torque 
with  01,  o2,  end  03,  we  obtain 


132  *  «*L 
Ip  Q  ’♦  PBI,  ■ 
Ip  B  -  ?QIg  •> 


Thus  we  can  write  the 

e 

P  - 


following 


MH 
MW 

three  equations  for  the  angular  acceleration: 


4  ■  ZH-PBIa 


If  ve  know  the  external  aoHsats  In  the  tea*  *  L,  «ZM,  and  XI,  we  can  solve 
for  P,  Q,  and  a,  by  integrating  tbs  respective  angular  accelerations. 

It  appears  that  ve  have  substituted  the  unknowns,  ML,  MM  and  JUT  for  tbs 
other  unknowns  P,  Q,  and  B,  hut  ve  eaa  determine  the  external  snsents  by  physical 
experiments,  whereas  It  nay  he  difficult  to  gat  the  angular  rates. 

We  now  have  six  aanodynaaic  quantities  which  aost  ha  evaluated.  Proa  the 
translational  equations  we  needed  three  serodyn— 1c  tsrws  D,  8,  and  L.  Bov  from 
the  rotational  equations  we  have  three  wore  aerodynamic  tarns  contained  in  J.L, 
ZM,  and  II.  Ve  will  discuss  these  shortly. 

Belatlve  and  Bxtarnal  Wad 


PToa  sane  source  (as  for  ample,  air- weather  at  K3MB)  we  need  to  obtain 
figures  for  the  external  wind.  Bw  values  which  ve  gat  for  the  external  wind 
will  ha  given  to  us  in  tarns  of  .directions  such  as  north  or  east,,  etc.  Proa  the 
given  direction  ve  aust  obtain  the  ocnponanta  along  each  of  our  ground  axes. 

Once  we  have  obtained  these  figures  in  the  ground  ayatsn  we  can  represent 
the  external  wind  as  in  the  following  equation: 

***  ■  Ifeil  ♦  V fe  ♦  ^«3 

fha  nsihod  for  dstamlhiag  the  relative  wind  eaa  be  eaaily  seen  from  the 
jfoUosring  analysis.  Suppose  ws  have  an  object  sitting  at  rest  on  the  ground 

TL 
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so  that  the  body-1  axis  is  parallel  to  the  ground-1  axis  and  the  external 
wind  is  bloving  parallel  to  the  ground^  axis.  Then,  as  fur  as  the  missile, 
itself,  la  concerned,  the  wind  acts  on  it  as  if  the  air  wars  still  aad  the 
missile  were  flying  backwards,  with  a  velocity  equal  in  Magnitude  to  that  of 
the  external  wind.  Shis  oan  be  sees  frea  the  following  diagram 


Thus  we  would  writs,  for  the  relative  wind  illustrated,  shove: 
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That  is,  If  we  know  the  direction  in  which  the  external  wind,  is  blowing,  then 
the  effect  on  the  missile  is  merely  the  negative  value  of  this  quantity. 

In  our  original  equation  for  the  external  wind,  ve  assumed  that  the  symbol* 
for  the  components  were  positive  nunbers.  The  sign  of  the  comments  is  deter¬ 
mined  by  the  direction  the  air  is  moving. 

The  first  thing  ve  must  do  vith  the  external  wind  vector  is  transform  it 
into  the  body  system.  That  is: 

■««  ■  "rfi  *  *  WjTj  -  VaJ?!  .  ♦  W«3t3 

The  components  ofthejsxternal  wind  in  the  body  system  are  found  by  taking  the 
dot  product  within,  and  respectively  and  utilising  the  transformation 
matrix.  This  gives: 

V^(,  rW4§.-S>W,&fc.*' Mfrh' 

We*,  z  Wx c*s4+- Wy  (-cos 4 sm 

Waa-  3,‘fci  +  W&'ij. + 
vVwls  sw'f’+VtyCes+GU't-  WgJiA+c*,  f* 

w«f|-  Wtt-%+ 

w4xj  3  - 


YGos&+*m4>svn6) 

+tfg(sir4imtM}*+c#+*i»+J 


The  velocity  of  the  missile  is: 

eeJb  ^-tr 

Va  •  ubx  ♦  vbg  ♦  Vbj 

noting,  that  the  relative  wind,  Wv,  is  the  velocity  of  the  minus  the  ex¬ 

ternal  wind,  we  have: 

Sr  •  hSL  ♦  ♦  vjb3 

where 

Vi  •  u  - 
«2  -  v  -  wex2 

Vv“  Vs 

and  we  have  completely  defined  the  relative  wind. 
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EXTERNAL  FORCES  AND  MOMENTS 


At  this  point  ve  have  derived  six  Independent  equations  describing  the 
missile's  motion: 


• 

<1  - 

*l/m  +  Rv  - 

s 

V- 

h/m  ♦  Pv  - 

Bu 

e 

TV*. 

Tl/a  ♦  Qa  - 

Pv 

e 

P  - 

n/in. 

Q  - 

(*M  -FRI,) 

1 

5. 

R  - 

UN  ♦  BJI,)  * 

The  mass,  a,  and  the  moments  of  Inertia  are  functions  or  curves  which  are 
obtained,  either  from  the  manufacturer  or  from  derivations  of  data  obtainable 
from  the  manufacturer.  Assuming  that  these  curves  are  available,  the  unknowns 
In  our  equations,  as  pet,  are  the  external  moments  for  the  rotational  equations 
and  the  external  forces  for  the  translational  equations.  The  external  forces 
arise  from  three  sources:  gravity,  aerodynamics,  and  thrust.  The  gravitational 
.and  thrust  forces  were  discussed  earlier,  so  the  aerodynamic  forces  are  the  only 
ones  which  we  need  to  discuss  In  greater  detail.  The  external  moments  are  of  two 
types:  aerodynamic  moments  and  moments  caused  by  thrust  misalignment. 

Earlier  we  defined  the  thrust  axis  system.  External  moments  will  arise  In 
the  body  system  due  to  the  position  and  direction  of  the  thrust  vector.  This 


There  ere  three  components  of  the  thrust  acting  In  the  body  system.  These  are 
found  by  transforming  the  thrust  vector  (as  was  dons  before)  Into  the  body  system. 
This  gives  us  the  equation: 

*  -  W  Vs*  ♦  *,  v 

Presuming  that  the  point  of  application  la  as  pictured  above,  the  vector  from 
the  origin  of  the  body  system  to  the  point  of  application  of  the  thrust  vector 
is: 


je 


The  vector  equation  for  the  aoasnt  due  to  thrust  ode &1  lament  la: 

Mw-txf, 

which  gives :  . 


Ve  can  write  the  three  terns  for  the  external  moment*  arising  from  thrust 
alsaligaaent  as: 

IAtH  r  LrT?/  ‘bt/b‘%,4‘Nr% 

Then  the  sum  of  the  external  moments  is : 

HLzlt+Lr 

where  the  subscript  "A"  iaplies  aerodynamic  meats. 

The  external  forces  can  be  written  as : 

Fi-Fn+T,  +  F,y> 

Fi  -  FiA  +  Ti  +  Fiir 
FtotTsf  Fjjr 

and,  from  both  sets  of  equation,  only  the  aerodynamic  functions  hare  not  been 
completely  defined. 

Bapresentatloa  of  Aerodynamics 

It  has  been  found,  after  many  years  of  experience,  that  the  aost  oomvUi 
aethod  for  representing  aerodynamic  farces  and  mo— arts  1st 

f>  =  $  P  */!sCf  *  )Scr  =-jrp  rt+scr 

\\-i9^5ZCtr  *slc*~  ftirfSicA 


wtosrs  Fp  is  an  arbitrary  aerodynamic  foroe,  1^  an  arbitrary 
ment,  P  tbs  density  of  the  air,  ^  the  air  prsesure,  6  and? 


surface  end  length,  respectively*  9  the  ratio  of  the  heat  capacities  of  the 
medium:  9  -  Cp/Cy  —  which  for  air  is  shout  l.hj  K  la  Mach  nunber,  Viy  the 

relative  wind,  Cp  and  CA  the  force  end  moment  coefficients,  end  q  .Is  the 
dynamic  pressure. 


Manor  times  in  the  past  (and  still  fairly  often)  we  were  given  the  aerody¬ 
namic  forces,  sad  a  position  on  thv  missile  called  center  of  pressure.  The 
center  of  pressure  is  tbs  point  at  which  the  aerodynamic  forces  act  on  the 
missile.  Since  the  missile  in  specs  Is  a  free- fro e  beam,  we  can  represent  the 
moments  and  forces,  independently,  as  acting  at  the  center  of  gravity.  The  sum 
of  the  external  forces  acts  through  the  c.g.,  end  the  sum  of  the  external  momenta 
act  as  torques  around  the  c.g.  fee  point  of  view  is  equivalent  to 
This  can  he  observed  from  the  following  diagram:  '  * 


The  equations,  as  defined  earlier,  assume  that  vc  bars  force?  acting  through 
the  center-of -gravity  and  that  we  have  moments  tending  to  rotate  tfc®  *uK3He  about 
the  c.g. 


The  air  density,  £  ,  in  the  force  and  moment  equations,!*  determined  frees  a 
"standard  atmosphere*.  V«  treat  the  points  on  the  curve  for  fee  we  did  the 
points  ou  the  curves  for  mass,  moment  of  inertia  and  thrust,  from  the  above 
aquations,  than,  only  Jfcch  number  and  the  coefficients  are  left  undefined.  Mach 
number  is  defined  as : 


where  V  is  the  local  speed  of  sound.  Ve  can  represented  as  a  function  of 
altitude.  Xt  la  generally  a  "standard"  curve,  m  is  the  air  density. 

Tbs  three  force  coefficients  are  defined  t&.tbe  wind  system  as  follows: 
Cp  -  dreg  force  coefficient 

Cg  -  aide  force  coefficient 

Cp  ■  lift  fore*  coefficient 


#  bene  rally,  8  is  the  area  of  the  ernes  section  of  tfca  missile  having  the  greatest 
Hamster  and  *  is  the  "moan  aercdynsmlr?  ohord." 
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or  if  they  are  defined  in  the  body  system  they  are  named : 

Cc  «  chord  force  coefficient 
Cy  ■  yaw  force  coefficient 
Cjf  ■  normal  force  coefficient 
The  moment  coefficients,  defined  in  the  body  system,  are: 

C£  -  rolling  moment  coefficient 
C  -  pitching  moment  coefficient 

a 

Cq  a  yawing  moment  coefficient 

Sot  ice  that  the  subscripts  are  capital  letters  for  forces  end  lover  case  letters 
tor  moments. 

It  has  been  found  in  most  cases  that  ve  can  represent  the  aerodynamics  as 
a  xfim  of  coefficients,  each  of  vhlch  is  a  function  of  two  variables,  one  of  thorn 
Bfach  number.  Instead  of  Mach  number,  we  sometimes  use  Reynold's  nuaber  for  low 
velocities,  end  ve  may  use  a  parameter  called  Tr oudo’s  number  for  very  high  al¬ 
titudes  where  the  air  cannot  be  considered  as  a  continuous  medium*  But,  as  far 
as  we  are  concerned,  our  independent  •jsaa&er"  will  be  Mach  mafcer.  for  most 
purposes,  lust  a.  fev  variables  a*  fuactioa*  of  Mach  masher,  will  suffice  to  de¬ 
scribe  the  aerodynamic  forces .  for  e»a$i£,  oC ,  the  eagle  of  attack  in  pitch, 
will  soeerttass  m$srm$sst  ivs  high  as  80$  of  the'  total  aerodynamic  pitch  force 
or  S  ,  the  sideslip  sx^ie  or  yawing  of  atta&k,  will  do  the  asms 

for  the  yaw  aemmi  or  jw  force ,  f&s  second  s#*t  important  parameter  (from  the 
standpoint  of  sisssiio  control  it  is  the  spat  important)  is  the  fin  deflection, 

S  .  addition  to  the  two  parameters  auctioned  above,  that  is  angles  of 
attack  and  fin  defbectisas,  ether  parameters  often  used  are  the  eguars  of  the 
angles  of  attack,  the  time  rate  of  ch sage  of  the  angle#  <sf  attack,  the.  piiciar*, 
yawing,  and  rolling  rates,  and  eo^irasixons  of,  say,  the  angles  of  attack  and 
fin  deflection  —  and  other  interdependent  parameters. 


Da  moet  representations  of  alee  lie  syat.as*  vs  'linearise*  the  aerodjnsmice* 
Ve  can' write  the  aerodynamic  coefficient,  say  C^,ae  a  function  of  its  various 

~  v-*j  {«* j  £ ‘ ' 

tha  derivative  of  this  function  and  assuming  independence  of  sash  vari¬ 
able,  wfe  retain: 


* 


In  standard  usage  the  notation  is: 


and  ve  can  vrite: 


If  vo  now  shift  the  burden  of  variation,  of  the  tern  due  to  change  in  Mach 
number  to  the  coefficients,  wt  can  write  the  above  as: 


Integrating  this  expression  gives: 

'  Cm-  Cm*<X+  Cnu  S f  «*+?»*  *+ C»*Q+  "  * 

which  gives  us  a  linearized  tenr  for  the  pitching  moment.  Although  the  above 
equation  looks  very  beautiful,  and  is  generally  quite  accurate,  it  is  not  exact. 
This  is  because  a  coefficient  is  not  necessarily  independent  of  the  separate 
variables .  To  the  contrary.  A  change  in,  say  will  perhaps  affect  Gm*  • 

That  is,  the  value  of  the  coefficient  is  not  the  same  for  different  values 

of  £  .  Therefore,  if  ve  want  greater  accuracy,  ve  would  have  to  include  equa¬ 
tions  of  tha  form:  etc.  As  was  said  before,  ve  are  generally  more 

accurate  if  ve  aesun  Independence  of  the  different  parameters .  than  the  actual 
accuracy  of  the  aerodynamic  coefficients,  theme e Ires .  The  above  parameters 
(which,  in  the  example,  were  shbwn  as  functions  of  the  pitching  moment)  are,  in 
general,  the  most  important  aerodynamic  parameters. 
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shows  a  change  in  <K  that  is  equal  to  Q.  _ 

With  the  advent  of  the  last  series  of  equations!  presuming  that  the  Man¬ 
ufacturer  or  designer  of  the  particular  vehicle  under  study  has  provided 
functions  for  the  aerodynamics!  then  11  the  unknowns; which  began  with  the 
original  force  equation,  are  now  too-  -  except  for  S  •  The  aerodynamic  func¬ 
tions  furnished  by  the  manufacturer  will  probably  not  be  in  exactly  the  fbna 
shown.  For  example,  seme  missiles  only  require  such  aerodynamic  functions 
are,  say,c(,  £,  and  Q.  Other  manufacturers  bring  in  completely  different 
parameters.  For  example,  the  aerodynamics  of  one  missile  use.  the  "bank  angle"  . 
nils  is  an  angle  which  lies  between  the  pitch  angle  of  attack  and  the  hypotenuse 
generated  by  the  side  slip  angle  and  the  pitch  angle  of  attack.  This  is  IUust 
t rated  in  the  following  diagram: 
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Notice  that  the  "hypotenuse"  is  actually  the  total  angle  of  attack*  OtT.  What* 
ever  form  the  aerodynamics  are  represented  in*  their  relationships  are  of  the  1 
form  shown  above. 

As  sons  tines  happens*  linearized  representation  of  the  aerodynamics  is 
not  sufficient.  (Repeatedly*  throughout  this  piper*  attention  has  been  called 
to  possible  sources  of  error*  approximations  used*  etc.  For  exsqple*  it  was 
pointed  out  that  the  surface  of  the  earth  is  not  actually  an  inertial  system 
nor  is  it  flat;  and  that  ve  had  to  find  a  set  of  a»s*  other  than  the  body 
axes*  to  correspond  exactly  to  the  principal  axes  of  inertia).  Generally* 
our  error  In  linearization  is  not  very  great.  If  it  is*  then  ve  oust  find 
other  methods  for  dete mining  the  aerodynamic* .  If,  for  example*  ve  represent 

as  a  curve*  and  if  it  Is  linearized*  it  will  be  a  plot  somewhat  as  follows 


M 

where  M  is  Mach  number.  If*  however*  were  also  a  function  of*  say*  0,  the 
roll  angle*  then  our  curve  would  have  to  be  "three  dimensional",  so  that  by 
moving  in  the  "third  dimension"  (corresponding  to  a  change  of  0)  we  would  find 
the  desired  coefficient  as  a  point  on  the  surface. 

TRAJECTORY  INFORMATION 

She  ultimate  purpose  of  an  analysis*  such  as  this*  is  to  determine  the  tra¬ 
jectory  of  the  missile  under  the  Influence  of  the  various  parameters.  As  was 
stated  earlier*  ve  made  our  velocity  vector  the  invariant  between  the  Inertial 
system  and  the  rotating*  or  body  system.  As  a  general  rule*  trajectory  Infor¬ 
mation  demands  the  range  vectors  as  a  function  of  time.  We  wish  to  know  the 
ranges*  but  ve  wish  to  know  these  ranges  Ir  the  ground  system.  Therefore*  ve 
will  Integrate  the  accelerations  in  the  body  system*  that  is*  u*  v,  and 
obtain  u*  v*  v*  the  velocities*  and  transform  them  into  the  ground  system. 

These  axe: 

V*  9i  +  YySi*’ Vsfj-  Vr  ubrhvbz  +  w 'tj 

Vx,  Vy-  and  are  found  by  forming  the  dot  product  vith  St*  3*.  *  *  respec¬ 

tively  and  utilizing  the  transformation  matrix.  To  obtain  the  three  ranges  in 
the  ground  system*  ve  now  Integrate,  these  three  components  of  velocity.  This 
gives  us  our  ground  range.  If  ve  keep  track  of  the  elapsed  time  during  com¬ 
putation*  ve  have  a  representation  of  the  ranges  in  the  ground  system  at  any 
instant  of  time. 
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FIN  DEFLECTIONS 


Assuming  that  the  analysis  Just  developed  is  for  a  missile  vith  guidance, 
then  there  is  some  form  of  external  "guidance  signal"  which  is  an  input  to  the 
guidance  circuitry.  Let  us  represent  the  "guidance  signal"  or  "guidance  cotasand" 
by  £c.  £e  may  be  a  radar  signal,  a  doppler  shift,  signals  from  a  stable  plat¬ 
form,  etc.  The  "fin  deflection"  may  not,  literally  be  an  angle,  but  may  be  a 
proportional  angle,  a  spoiler  effect,  a  dwell- time,  etc.,  but  in  any  case  we  can 
represent  the  relationship  of  the  "delta  cosmand,"  £c,  to  the  "delta  achieved", 
by  means  of  a  transfer  function.  Thus : 


cSc 

KG(s) 

/n/s  SIL£ 

^QNS, 

r 

Vith  cJc,  as  yet,  our  only  unknown,  we  can  represent  the  entire  missile  airframe 
(with  the  guidance  represented  as  S'c.  XB(s )  «  £ )  by  means  of  block  diagrams . 

These  block  diagrams  help  us  bo  visualize  the  interrelationship®  and  interdepen¬ 
dences  of  the  equations  which  we  have  been  deriving.  On  the  following  page  is 
a  block  diagram  of  our  system  of  equations. 

The  mathematics  for  every  missile  system  will  be  basically  similar  to  the 
.  mathematics  we  have  derived.  This  system,  described,  is  far  from  being  complete. 
As  yet  we  have  not  mentioned  guidance,  which  will  involve  such  things  as  seekers, 
targets,  radars,  preprograms,  etc.  What  we  have  done,  however,  is  mathematically 
describe  the  missile  dynamics,  and  indicate  methods  for  more  exact  analysis  if 
needed. 

TRANSFER  FUNCTIONS 

In  the  previous  analysis  we  had  but  one  unknown,  the  fin  deflection,  c T  . 

The  "fin  deflection"  represents  either  the  actual  fin  position  as  an  angle,  or 
a  "dummy"  parameter  representative  of  the  fin  deflection.  If  etc  is  bhe  command 
to  the  missile,  we  can  then  represent  the  process  between  <fc  and  (Tut  trans¬ 
fer  function,  KG(s).  A  transfer  function  is  merely  a  Laplace  transform  of  the 
hardware  in  the  missile  with  zero  initial  conditions.  Transfer  functions,  as 
such,  will  not  be  taken  up  in  this  paper.  Suffice  to  say,  there  exist  a  number 
of  good  texts  on  the  Laplace  transform  and  transfer  functions*. 


*  For  example,  see:  Thaler.  0.  S.  and  Brown,  R.  0.,  Servomechanism  Analysis, 
McGraw-Hill,  1953* 
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For  any  type  of  problem  la  which  a  texgetia  used,  the  target  poaitiaa  and 
velocity  will  be  represented  in  ground  coordinates.  She  following  diagram 
Illustrates  the  target  range  vector  and  its  caqponesrU  along  the  three  ground 
arts  j  ■ 


Thus,  the  target  range  vector  Is: 

*^T  “  +  %y92  * 

If  we  alio  know,  the  target  velocity  we  can  represent  it  as: 

She  coordinate  system  fixed  to  the  surface  of  the  earth  (whether  it  is  assumed 
to  be  the  inertial  system  or  whether  it  is  rotating  and  the  center  of  the  earth 
is  taken  as  the  inertial  system)  represents  a  convenient  way  (in  fact,  probably 
the  only  way)  which  ve  can  introduce  target  ranges  and  velocities  late  our 
system  of  equations. 

At  'tlme*  ithe*  analyst  is  allowed  to  bs  fairly  arbitrary  la  his  choice  of 
ranges  and  velocities  of  the  target  ••  insofar  as  these  choices  represent  values 
whlohu.can  be  realised.  Many  times  It  is  convenient  to  have  a  target  with  * 
velocity  along  one  of  the  ground  axes  but  no  velocity  component  any  other 
axis .  This  eases  the  complexity  of  the  problem.  When  this  situation  arises 
(that  is,  when  a  target  has  a  constant  velocity,  perhaps,  or  even  a  cHfttifrt 
range)  it  Is  often  desirable  to  orient  the  ground  system  so  that  the  target  " 
flies"  parallel  to  'one  of  the  unit  vectors. 
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the  sotted  for  obtaining  the  individual  components  of  the  target  range 
and  velocity  will  not  be  the  sane  ae  ve  have  discussed  before.  Previously  ve 
transformed  from  one  coordinate  system  to  another  by  means  of  an  Euler  angle 
sequence.  This  method  vas  used  because*',  we  wished  to  rotate  one  set  of  axes 
(in  this  case,  the  ground  axes)  so  that  their  final  position  would  be  aligned 
with  another  set  of  axes  (in  this  case,  the  body  axes).  Ve  can  assume  that  the 
target  is  a  point  in  space.  Therefore,  ve  have  a  vector  which,  if  it  is  the 
range  vector,  reaches  from  the  origin  of  a  ground  system  to  a  point  which  ve 
call  the  target.  This  vector  is  then  represented  as  having  three  components 
parallel  to  the  ground  axis  system.  In  order  to  determine  the  components  of  the 
•  target  range  (or  velocity)  we  will  use  a  method  called  "rectangular  projection". 
First  we'  obtain  a  projection  of  the  vector  in  the  3f,  -  91  plane,  an  intermediate 
range  vector,  ,  as  in  the  above  diagram.  This  length  is  found  by  rotating 
thro  ja.  the  target  elevation  angle,  ftp,  as  in  the  following  diagram: 


Rrz 


which  gives  for  Rrpz  and  the  values : 

\ 

ft-rj  =  J  ftTl  cos&f 


RtzSI^tI  5,n'&r 

The  following  diagram  illustrates  this  new  component  of  the  original  vector; 


If  ve  Introduce*  target  azimuth  angle,  yJr  ,  then  the  other  two  components  of 
the  target,  la  the  -  gg  plane  -re: 

Rr  y  =  /  Rr  I  COS  0r  CM  %. 

RTa  = j  COS  O'  sinipr 
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and  ve  have  completely  represented  the  target  position  vector* 

The  projection,  as  defined  above,  assumes  that  our  intermediate  vector 
is  projected  onto  the  horizontal  plane,  that  is,  the  ^  plane.  One 
may  vender  if  it  wouldn't  be  as  correct  to  project  the  original  range 
vector  onto  a  vertical  plane  such  as  the  ^  -  g?  or  gg  -  plane,  and, 
from  the  intermediate  vector  so  projected,  get  the  components  of  the  range 
vector.  Mathematically  and  geometrically,  either  method  is  equivalent. 

The'  first  method  corresponds  to  a  0  projection  followed  by  a  Y  projection* 
the  second  to  a  'Y,  9  projection. 

The  limitations  imposed  on  the  choice  of  the  projection  depends  on 
which  angles  are  measured.  In  the  above  Instance,  and  Qj,  are  angles 
of  a  type  which  a  radar  system  might  measure.  The  vertical  ingle,  <5j,  is 
measured  by  a  radar  from  the  horizontal  plane,  that  is,  tha1^  - *gg  plane, 
and  the  horizontal  angle,  'f’m,  is  measured  in  this  same  horizontal  plane. 

In  such  a  situation,  ve  ate  forced  to  use  a  vertical  projection  first* 

This  can  be  demonstrated  by  reference  to  the  following  figure  i 


If  ve  project  flret  onto  the  X-Z  plane  through  the  angle  then  for  the 
azimuth  angle  ve  have: 

If  ve  projeot  first  onto  the  X-Y  plana,  then  the  sine  of  the  aslauth  angle, 
Tg1** 


furthermore,  projecting  first  on  the  X-Z  plane  gives 


for  the  elevation  angle,  and 


if  we  first  project  the  range  vector  onto  the  horizontal  plane. 

Obviously  our  projection  method  depends  upon  vhlch  angles  were 
measured.  Radar  angles  necessarily  proceed  through  a  0  projection 
followed  by  a  ‘f . 

At  this  point  it  is  difficult  to  say  what  will  be  done  with  the  tar* 
get  range  and  the  target  velocities  once  they  are  expressed  in  the  ground 
coordinate  system.  Each  kind  of  missile  has  a  different  guidance  philosophy 
and  will,  therefor*,  utilise  these  parameters  in  different  ways.  But,  in 
every  case,  the  guidance  philosophy  must  use  them  in  one  way  or  another. 

We  will  not  attempt  to  discuss  guidance  now,  but  will  Illustrate  bobs  types 
later.  Suffice  to  say,  we  can  transform  the  target  velocity  vector  Into  > 
or  missile  body  coordinates  if  we  wish;  wo  can  use  the  target  range  for  a 
radar  ground  station;  or  vs  can  utilize  both  these  vectors  in  several  dif¬ 
ferent  ways.  If,  pbrhape,  the  guidance  of  the  missile  needs  the  distance 
of  the  missile  to  the  target,  then  vs  can  find  the  relative  distance  between, 
target  and  missile  simply  by  subtracting  the  components  in  the  ground  ayw  , 
and  then,  if  we  wish  to  express  these  relative  components  In  missile  body 
coordinates,  we  can  transform  them  inaed  lately  with  the  aid  of  the  transfor¬ 
mation  matrix. 


As  a  closing  paragraph  to  this  subject  we  might  mention  two  situation* 
In  which  we  would  need  to  transform  these  vectors.  For  example.  If  the 
■  guidance  of  the  missile  depends  upon  the  intensity  of  seme  source  of  radia¬ 
tion  relative  to  a  guidance  head  in  the  missile,  such  as  infra-red  radiation, 
then  the  brightness  of  the  radiation  at  the  surface  of  the  detector  would 
depend  upon  the  distance  of  the  missile  from  the  source.  In  this  case  we 
would  like  to  know  how  far  the  missile  is  from  the  source.  After  finding 
this  distance,  we  would  need  to  transform  it  into  missile  coordinates  so 
that  we  could  find  the  relationship  between  the  direction  of  the  detector 
and  the  missile  heading  direction  —  and  determine  tracking  error,  capability, 
etc.  As  another  example,  the  maneuvers  of  the  missile  might  depend  upon  the 
relative  velocity  of  the  target  with  respect  to  the  missile.  In  this  case 
we  would  find  the  relative  velocity  between  the  target  and  missile  and  trans¬ 
form  the  resultant  components  into  the  body  system.  Ibis  would  give  desired 
information  to  the  missile. 

GYROS 


Before  beginning  the  discussion  of  gyros  vs  would  like  to  Insert  at 
this  point  a  demonstration  of  a  statement  which  we  will  need  later  in  the 
discussion.  We  recall  that  earlier  vs  derived  a  transformation  matrix 


,  between  the  ground  and  body  systems .  This  particular  matrix  was  called  an 
"orthogonal"  matrix.  One  property  of  this  matrix,  vhlch  is  Jjgortant  in 
the  following  discussion,  is  that  each  element  of  tbs  matrix  is  equal  to 
its  co-factor.  For  example,  a  matrix  using  the  y,  i,  and  0,  In :  . 


asanmMfrftV’fl* am# cai #**>>' 
•  ots+s*?  casfeesf' 

as9p*tsmf\ -**m#c*  ?+&***? 


«*♦  I 

I 


% 

b. 

s. 
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Bier  ten  In  the  upper  left-hand  comer  of  the  above  matrix  has,  ae  ite  co- factor, 
the  detemlnant  found  by  deleting  the  flret  row  and  first  colum.  If  we  solve 
this  detemlnant  we  get  the  following  tern: 

cosjcosys  zcostyc&eccsi'+jirfyeasws? 

which  Is  the  first  tern  In  the  Matrix.  This  sane  type  of  operation  can  be  shown 
for  every  element  of  the  Matrix. 

In  a  Missile  which  obtains  its  position  by  naans  of  gyros,  generally’  two 
tvo-degree-of- freedom  gyros  am  employed.  These  gyros  am  "slaved”  In  order  to 
measure  all  of  the  angles.  The  foUovlxg  two  diagrams  Illustrate,  this.  Bach 
diagram  shove 'a  schematic  of  a  gyroscope  having  a  pitch  outside  global. 


Gyroscope  a,  cf  course*  may  have  fixed  outside  gimbals  other  then  pitch. 
There  are  six  possibilities  in  all:  two  pitch,  two  yaw,  and  two  roll.  The 
second  angle  further  differentiates  between  types  of  gyros,  there  being  two 
possibilities  for  each  type,  as  is  Illustrated  In  the  abore  diagnose. 

At  the  center  of  each  gyroscope  is  a  spinning  mass  whose  axes  of  rotation 
remain  fixed  In  space.  If  the  initial  direction  of  these  axes  is  known  with 
respect  to  the  ground  coordinates,  then,  within  the  missile  Itself,  is  a  rector 
which  retains  a  fixed  direction  is  ground  coordinates.  The  fixed  axis  is  spin¬ 
ning  very  rapidly,  but  by  integration  of  the  spinner  angular  velocity  an  angle 
about  this  axis  can  be  determined.  We  will  not  actually  use  this  angle  because 
it  wr  d  not  be  accurate*  We  will  use  only  pick -off  angles  of  the  two  outer 
glmba.  rings.  We  can  see  that  a  sequence  of  angles  from  the  Intenaost  spinning 
mass  to  ths  outer  glmbdl  ring,  which  is  fixed  to  the  missile,  will  give  a  se¬ 
quence  of  angles  Which  determine  the  missile  position  in  space  with  respect  to 
the  ground  coordinate  system. 

We  cannot  "arbitrarily"  choose  the  Euler  angle  sequence  of  the  angles  in 
the  gyrosoope  as  we  did  in  the  transformation  between  the  ground  and  the  body. 
Once  the  gyroscope  is  mounted,  the  sequence  of  angles  is  fixed,  The  necessary 
Euler  angle  sequence  can,  however,  be  determined  by  inspection. 

The  pick-off  angles  of  the  gyroscope  uniquely  determine  the  orientation 
of  the  missile*  We  can  demonstrate  this  very  easily  In  the  following  aannar. 

The  angular  velocity  vector  can  be  represented  as: 

00  a  Pb\t  Qbg  t  fits  s  SfSz  *  +  Vji?j 

where  the  subscript  "g"  represents  the  gyroscopic  pick -off  angles.  The  gyro¬ 
scopic  pick-off  angles  are  the  integrals  of  the  following  quantities: 

£>3<sL  -  -  - 

Yj-  fS*S  +Q%-Z+Rb3''?j 

The  integrals  of  the  equations  show  that  the  pick -off  angles  are  unique  since 
they  are  determined  entirely  by  the  orientation  of  the  missile  in  j-race. 

In  diagram  "I",  the  sequence  (as  can  be  seen  by  inspection)  is  a  yaw  atgl= 
(with  reference  to  the  spinner)  then  a  roll  angle  followed  by  a  pitch  angle.  In 
diagram  "II",  the  eequence  is  a  roll  angle  (corresponding  to  the  spinner) 
by  a  yaw  angle  and  than  by  a  pitch  angle.  These  two  sequences  are  from  the  "in¬ 
side  out”.  We  can  write  the  two  transformation  equations  in  matrix  notation  ast 

(b;)s  MMMW* M(J£)ST 

Wmn  the  M*e  represent  transformation  matrices  in  the  indicated  modes. 
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Let  us  determine  the  transformation  for  the  two  types  of  gyros  illus¬ 
trated.  For  diagram  Bl",  ve  will  indicate  the  angles  by  writing  subscript 
"V”  for  each  angle,  m  the  second  one  we  will  use  a  subscript  "H”  for  es^h 
angle.  These  subscripts  mean  that  the  first  gyroscope  is  "vertically  fixed" 
and  the  second  is  "horizontally  fixed".  That  is,  the  spinner  position  is 
vertical  or  horizontal,  respectively.  Fran  the  vantage  point  of  our  previews 
knowledge,  we  can  choose  the  intermediate  matrices  for  pitch,  yaw,  and  roll. 
We  will  indicate  these  matrices  below: 
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Having  these  three  matrices  available  and  choosing  the  subscript  "v"  for  th* 
first  transformation,  and  the  subscript  "H"  for  the  seconi  transformation  vs 
obtain: 
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Now  that  we  have  the  actual  values  of  the  tnmaformtion  mtriees,  we  will 
rewrite  the  last  two  transformation  matrices  in  a  slightly  easier  to  handle 
notation.  Let  the  following  two  matrices  represent  the  "vertical'"  and 
"horizontal"  matrices  respectively: 


The  spinner  has  special  significance.  The  orientation  of  the  spinner  axis 
remains  fixed  throughout  the  time  that  the  gyroscope  is  in  operation.  (This  is 
not  quite  true  since  we  are  as  summing  that  the  earth  is  flat,  and  the  orienta¬ 
tion  of  the  gyro  will  remain  fixed  with  respect  to  space  rather  than  the  earth's 
surface.  It  will  appear  that  the  direction  of  the  vectors  changes  by  an  angle 
equal  to  the  change  of  the  normal  vector  of  the  earth's  surface  due  to  Its  cur¬ 
vature.  However,  for  this  discussion  we  will  assume  this  angle  is  small  enough 
to  ignore).  If  the  missile  is  originally  "on  the  ground"  and  its  vectors,  the 
body  axes,  are  aligned  parallel  to  the  ground  axes,  we  can  make  the  following 
statement: 


That  is,  the  third  axis  of  the  vertically  fixed  gyroscope  always  points  in  the 
direction  of  the  third  ground  axis,  and  the  third  axis  of  the  horizontally  fixed 
gyroscope  always  points  in  the  same  direction  as  the  first  ground  axis.  Let  us 
now  write  the  transformation  matrix  which  we  originally  had  between  the  ground 
and  body  system  in  short  hand  notation  as  follows : 


/  *H  *12  *13  V| 
*21  *22  *23 


\*3i  *32  *33/ 
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Then  we  can  write,  after  replacing  the  vector  by  g.  in  the  first  transforma¬ 
tion  -matrix  for  the  gyroscopes :  * 

T, .  q  -  V 

H.=v„ 

-'v« 

and  from  the  transformation  matrix  between  body  aad  ground  we  obtains 

lr, "£?j  r  dn 

V^j =a»» 

V  ?|  "  ^-Sl 

•o  we  can  write: 

V(3  =  A,j 
Vjj  r  dii 

Yst =  a»J 

la  the  traaafozwMitlon  matrix  for  the  ”11”  gyroscope,  we  can  replace**^  by  g^. 
This  gives  us  the  three  equations: 

Vt(=H„ 

=H»i 

and  from  the  treos formation  matrix  between  body  end  ground  we  obtains 

V?, =  a,( 

elk  -4 

U#*9fj  -  du 

tj-f,  =  a„ 


» 


Therefore,  we  can  write: 


H„  = 

H3i  - 

H  is  s  &3I 


and  .re  now  write  for  the  transformation  matrix  between  body  and  ground: 


Utilizing  the  properties  of  orthogonal  matrices  explained  at  the  beginning 
of  this  section,  we  can  solve  for  the  middle  column  of  the  above  matrix. 
This  is: 


atl  * 


H2 1  v„ 

H31  Vj* 


H 
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We  have  now  completely  determined  the  angles  between  the  body  and  the  ground 
frame  in  terms  of  the  pick-off  angles  of  the  gyros. 


In  the  transformations  between  the  gyro  system  and  the  body  system,  the 
vertically  fixed  gyro  used  the  sequence:  yaw,  roll,  and  pitch.  The  horizontally 
fixed  gyro  used  the  sequence  roll,  yaw,  and  pitch.  It  was  pointed  out,  when  we 
were  performing  the  transformation  between  body  and  ground,  that  the  second 
angle  must  never  approach  very  close  to  90°.  This  was  because,  whan  determin¬ 
ing  the  roll rates  of  the  missile,  it  was  necessary  to  use  the  secant  of  the 
second  angle.  If  the  second  angle  ever  approached  9 0°,  then  the  value  of  the 
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entire  tern  would  become  ir-lofonnijoaiit.  An  analogous  situation  exists  with  the 
gyro  angles.  If  the  :.*coi;l  ancle  goes  to  90°  then  a  situation  arises  calle'’ 
"gimbal  lock".  Therefore,  In  the  vertical-fixed  gyro,  the  roll  must  never 
approach  90°,  and  in  the  horizontal  -  fixed  gyro,  the  yaw  must  never  approach 
90°.  If  they  do,  and  gimboi  lock  occurs,  then  it  is  no  longer  possible  to 
know  the  orientation  of  the  missile.  The  physical  reaction  of  the.  gyros  is 
often  a  complete  l0O°  turn  of  the  entire  gyro  setup  —  and  e  new  set  of  orienta¬ 
tion  vectors. 

With  our  previous  knowledge  of  transformation  between  systems,  we  were  able 
to  obtain  the  complex  transformation  equations  for  utilization  of  the  gyro  pick- 
off  angles  very  readily.  If  we  tried  to  visualize  the  complicated  interrela¬ 
tionships  of  the  pick-off  angles  to  each  other  and  also  the  interrelationships 
to  the  transformation  angles  between  the  ground  and  the  body,  it  would  have  been 
impossible  to  solve  the  equations.  But  with  the  utilization  of  coordinate  . 
systems  and  transformation  matrices  between  these  systems,  we  can  find  the  values 
in  a  very  simple  manner.  It  is  not  until  we  try  to  evaluate  the  terms  derived, 
such  as  the  V's,  the  H’s,  and  the  a's  of  the  matrices  that  the  equations  begin 
to  look  complicated.  Even  then  they  are  not  difficult,  they  merely  have  a  great 
many  terms. 

A  separate  section  will  not  be  introduced  for  an  explanation  of  stable  plat¬ 
forms.  The  main  difference  between  the  representation  of  stable  platforms  sal 
gyroscopes  is  that  a  stable  platform  (which  is  stabilized  by  means  of  gyroscope^, 
has  available  three  pick-off  angles.  In  our  analysis,  using  gyroscopes,  we  re¬ 
quired  two  of  them  in  order  to  determine  the  three  angles.  Bad  we  decided  to 
pick-off  the  angle  for  the  spinner  and  evaluate  its  size  (by  integration  of  the 
'  spinner  angle),  we  could  have  done  a  very  similar  thing  as  a  stable  platform 
would  do.  As  a  matter  of  fact  we  "pretended"  in  our  transformation  between  the 
gyro  system  and  the  body  system  that  we  actually  had  this  fixed  direction  angle. 
We  did  not  use  this  fixed  angle,  but  determined  the  orientation  angles  of  the 
body  with  respect  to  the  ground  from  only  the  outer  two  pick-off  angles.  If 
anyone  who  is  using  these  notes  finis  it  necessary  to  determine  the  equations 
relating  a  missile  system  to  a  stable  platform,  then  he  need  only  remesfcer  that 
he  transforms  from  the  stable  platform  to  the  body  system  through  three  anglec . 
And  also  therm  is  no  second  gyroscope  to  worry  about.  Therefore,  the  stable 
platform  is  more  simple,  in  essence,  than  the  gyroscopes  which  ws  Just  inves¬ 
tigated. 

.GVIPAO 

Included  in  the  acccepenylng  sections  are  typical  examples  of  guidance  and 
some  of  the  nomenclature  one  is  likely  to  encounter. 


GUIDANCE  AND  CONTROL  SYSTfcM  FOR  A  TYPICAL  RADIO  CGMIARD  MISSIIE  SYSTEM 

By:  Barry  McJunkin 

In  such  a  system,  the  chief  inputs  Into  the  ground  computer  are  the 
position  coordinates  of  the  target  and  of  the  missile.  These  coordinates, 
for  the  typical  case,  are  range,  elevation  angle  and  azimuth  angle  as 
measured  by  radar  instrumentation. 

The  ground  computer  transforms  these  coordinates  into  rectangular  co¬ 
ordinates  in  a  common  frame.  The  target  coordinates  are  time  differenti¬ 
ated  to  determine  present  target  velocity,  then  the  time  of  missile  flight 
to  a  predicted  intercept  of  the  target  is  determined  from  use  of  present 
target  position,  present  target  velocity,  and  stored  missile  capability 
curves.  Of  course  the  computations  are  modified  to  take  care  of  the  offset 
r  ti  ns  lent  time  needed  for  most  effective  use  of  the  warhead. 

Other  functions  necessary  for  efficient  steering  of  the  missile  are 
generated  from  their  established  dependence  on  the  coordinates  of  the  pre¬ 
dicted  burst  point.  These  functions  include  such  items  as  initial  dive 
order,  glide  bias,  and  missile  velocity  correction.  The  dive  order  is 
composed  of  acceleration  orders  to  the  missile  causing  it  to  turn  from  Its 
initial  nearly  vertical  trajectory  (given  during  the  boost  phase  which  is 
not  controlled)  into  a  senrl-parabollc  path  towards  the  predicted  Intercept 
point.  The  glide  bias  orders  are  commands,  ceparate  from  steering  orders, 
causing  the  missile  to  always  have  a  certain  fraction  of  1  g  acceleration  in 
a  nearly  vertical  direction.  The  missile  velocity  correction  added  to 
"present"  missile  velocity  during  steering  provides  more  accurate  determina¬ 
tion  of  tlme-to-go  and  consequently  of  the  steering  orders. 

The  steering  orders  are  determined,  essentially,  as  error  quantities 
proportional  to  the  differences  between  the  required  closing  velocity  coor¬ 
dinates,  and  th.  actual  closing  velocity  coordinates.  The  required  closing 
velocity  coordinates  are  the  differences  between  missile  position  coordi¬ 
nates  and  corresponding  target  position  coordinates,  divided  by  the  time-of- 
flight  to  the  predicted  intercept.  The  "present"  target  altitude  is  modified 
to  furnish  an  aiming  point  above  the  target  for  most  of  the  trajectory,  and 
thus  provide  greater  missile  end-game  maneuverability. 

During  a  large  part  of  the  trajectory,  the  remaining  time-of -flight  is 
continually  recalculated  in  terms  of  present  missile  and  target  velocities. 
This  function,  called  tlme-to-go,  replaces  the  function  initially  termed 
time-of-flight.  The  ve  locities  of  missile  and  target,  as  well  as  target 
accelerations,  are  determined  by  analog  type  differentiation  devices.  The 
actual  closing  velocity  coordinates  are  found  by  •  imply  subtracting  the 
iargei  velocity  coordinate  from  corresponding  missile  velocity  coordinates. 

Of  course  these  closing  velocity  coordinates  are  modified  by  the  missile 
velocity  correction  function,  and  by  present  target  acceleration  components. 

The  error  quantities,  described  above,  are  basically  velocity  differ* 
enoee  in  the  fixed  rectangular  ground  frame  associated  with  either  the  launch¬ 
ing  point,  the  missile  tracking  radar,  or  the  target  tracking  radar.  These 
quantities  art  converted  into  components  in  a  rectangular  frame  nearly  paral¬ 
lel  to  the  "pseudo"  missile  frame,  —  a  fixed  frame  attached  to  the  missile. 
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This  conversion  is  accompli* had  by  a  yaw,  pitch,  yaw  sequence  «s 
follows: 


First  rotate  about  gn  through  angle  A,  where  A  is  the  azinuth  angle  of 
corouted  intercept  point. 


Next  rotate  about  "g^  through  a  "clirib  angle"  9j(Q  will  be  derived  later) t 

& 

a  *er 

Z,x  _k 

‘L=^c<,s*,'%s/n* 

*»n*  43i,C9f9 


s»  ^ 

Finally  rotate  about  g^  through  a  "turn  angle"  (if  will  be  derived  later): 

6;  =  #,%cajV 5m  y 

^  +3ifcc«v' 


4* 


Jcsfelnlng  these  in  matrix  font  we  get  the  transformation  matrix  below  for 
the  peeudo-alsslla  frame:  .  .  « 

I  o\  I C OS  fCCSft - $\H+tCi*Si»A 

hi  I  =  [  -SififQ&A  -  COJ  rc*S4Siflti  ~£,nVMA+Cv+oai+c6sA  Mt/*0  I 
\T>lj  \  Sin9SimA  -jtntoesA  <*** /  \  1 
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where  V  ■ 


|VmJ  -  | V  +  vy2  4  Vs2 
taking  the  dot  product 2  with  g} ,  g^  end  gy  respectively,  we  obtain: 

Vbg-%  =  Vx  =  V(-SlNfCOSH -cos?  cose  SINM 
Vbr$  =  v/--  V(-SIVfS/A/H  +  cosy  cose  COSfl) 
ViW  =  V*  -  VCOSfS/N©. 

If  we  multiply  (l)  by  cos  A  and  (2)  by  sin  A  and  add  ve  obtain: 

VxC05Ai-V/5lNA  =-V5m 

or  ,ln»=-(Vx  COSfl 

and  If  we  multiply  (l)  by  sin  A  and  (2)  by  cos  A  end  subtract  (l)  from  (2) 
then: 

w  -Vx$iNft  i-tyCo/)  ^vcosycoso- 

Multiplying  (k)  by  sin  0  and  (3)  by  cos  0  and  subtracting  (4)  from  (3)  ve 
get: 

V/i iSlQSINe  -Vy  COSH  SIN  9  +Vj!  coso»o 

and  dividing  by  cos  0: 

tan  0  (Vx  sin  A  -  Vy  cos  A)  ♦  V  -  0 

or  ““  ®  ■  Vx  SllOT-  Vy  COSfl 

Since  the  missile  Is  stabilized  at  a  known  angle  of  roll  about  lie  longitudi¬ 
nal  axis;  and  with  respect  to  the  plane  containing  the  velocity  vector,  the 
conversion  of  error  quantities  into  the  "pseudo"  missile  frame  can  be  ccaf- 
pie  ted.  The  lateral  control  directions  of  the  missile  nearly  coincide  with 
the  lateral  axes  of  the  pseudo  missile  frame. 

Tiae-to-go  is  continually  determined  during  a  large  part  of  flight  as 
previously  mentioned,  nils  Is  done  in  such  a  way  that  only  lateral  ocemaade 
to  the  missile  are  different  from  zero,  since  normally  there  is  negligible 
control  in  the  longitudinal  direction  through  use  of  eleven  oootrol. 

Finally  the  lateral  error  quantities,  which  have  dimensions  of  velocity, 
are  converted  to  dimensions  of  acceleration  t) trough  division  by  tiae-to-go; 
then  these  are  scaled  to  the  proper  order  magnitudes  through  multiplication 
by  a  constant  factor.  The  orders  are  next  sent  to  command  circuits  which 
convert  the  orders  to  coded  electrical  pulse  groups.  These  groups  may  be 
generated  intermittently  and  alternatively  at  audio  frequency.  The  radar 
transmitter  accepts  these  pulse**  orders  and  transmits  the  information  via 
an  electromagnetic  link  to  the  mis •  lie  receiver. 
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'live  miss  lie  receiver  recognizes  properly  coded  signals,  reconverts  these 
to  pulse  groups,  and  seeds  these  to  a  signal  data  converter.  Each  time  a 
pulsed  group  is  recognized,  a  hold-off  signal  may  he  issued  to  burst  control 
circuits  and  an  activating  signal  sent  to  the  missile  beacon.  If  insuffi¬ 
cient  autopilot  control  signals  axe  received  the  hold-off  pulse  may  be  dis¬ 
continued  and  the  detonation  circuits  become  operative. 

Km  signal  data  converter  usually  recovers  In  dc  analog  form  the  ground 
computer  acceleration  orders  to  the  missile.  These  orders  are  sent  to  the 
autopilot  control  systems.  The  lateral  acceleration  control  systems  are 
usually  identical  in  form,  each  exercising  control  In  a  plane  perpendicular 
to  the  hinge  line  of  the  control  fins  of  the  system.  For  a  supersonic 
missile  the  lateral  control  direction  Is  nearly  normal  to  the  longitudinal 
aria  of  the  missile. 

A  lateral  control  system  may  be  comprised  of  a  steering  amplifier  with 
associated  Input  mixing  networks,  a  hydraulic  system  to  control  servo  power 
sa  actuator  assembly  driven  by  this  power,  a  mechanical  linkage  to  the  con¬ 
trol  elevons,  and  several  feedback  signals  to  the  steering  amplifier.  These 
feedbacks  may  Include  signals  representing  the  missile  acceleration  achieved 
in  the  desired  direction,  missile  angular  velocity  in  the  control  plane,  and 
average  control  fin  position.  The  remamnd  signals  are  compared  to  these 
feedback  signals  which  are  developed  by  missile  flight  instruments.  The  net¬ 
works  associated  with  the  steering  aapllfler,  and  into  which  the  orders  and 
feedbacks  are  Inputs,  compare  somewhat  continuously  the  present  response 
conditions  of  the  missile  to  the  present  steering  command.  The  result  Is  a 
smooth  response  to  orders,  and  a  preservation  of  missile  stability  throt^h- 
oot  the  flight. 

The  roll  control  system  operates  very  much  like  one  of  the  lateral 
control  systems,  excepting  that  the  roll  control  signal  arises  from  the 
missile  rolling  more  or  less  than  the  desired  roll  stabilization  angle,  as 
described  above.  Useful  feedback  signals  may  be  angular  roll  velocity  and 
roll  eleven  position. 

Efficient  mi se lie  systems  any  utilise  lateral  control  elevons  to  achieve 
roll  control  as  well  as  lateral  acceleration  control.  For  such  systems  the 
mechanical  linkages  are  arranged  to  that  average  control  eleven  position  is 
correct  for  the  deeirtsd  lateral  acceleration,  while  the  differential  between 
opposite  fin  deflections  is  oorrect  for  achieving  roll  control. 

The  following  diagrams  illustrate  a  typical  radar-command  roll-stabilized 
system  ae  described  above. 
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The  folioving  diagram*  illustrate  a  typical  radar  ccasnand  roll 
stabilized  system  as  above  described: 
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SEMI -ACTIVE  HOMO©  GUIDANCE  SYSTEM 


tfyacyragi rerw&gpw^  »•>  a-**-*  •**  •«  ■  *" 


!Rr  Merle  Lesoing 


In  »  setr-i -active  basing  guidance  system,  &  proposed  target  is  radiated 
with  RF  (radio  frequency)  energy  from  setae  source  external  to  the  missile,  and 
the  missile  utilizes  the  energy  reflected  off  of  the  target  to  derive  guidan-.'^ 
information,  For  example,  a  ground  “based  Ctf  ( continuous  wave)  illuminator  track? 
the  proposed  target  and  radiates  it  with  RF  energy*  The  missile  receives  the 
;-ef.V  '“-ad  energy  through  a  gimballed,  forward -looking  antenna  dish.  The  antenna 
has  i:  .rutating  beam  which  allows  the  missile  to  derive  tracking  errors  which 
cause  the  dish  to  always  point  toward  the  target. 


|  In  order  to  detect  and  to  continue  tracking  one  specific  target  out  of 

I  many  possible  targets  that  may  “be  within  the  “beam  width  of  the  Illuminator,  the 

I  internal  circuitry  may  “be  equipped  with  a  very  narrow  velocity  "gate"  which  it 

|  uses  to  resolve  a  single  target.  Thus  the.  antenna  tracking  errors  can  them  be 

l  derived  for  this  one  target.  The  velocity  gate  makes  use  of  the  doppler  shift, 

!  which  is  the  shift  in  frequency  of  a  signal  due  to  the  radial  velocity  between 
the  transmitter  and  the  reflector,  A  "rearward-looking”  entemia  in  the  missile 
heterodynes  a  sample  of  the  transmitted  energy  from  the  Illuaitntcr  with  the  re¬ 
flected  target  signal  received  at  the  front  antenna*  The  difference  in  fre¬ 
quency,  called  the  "doppler,  is  proportional  to  the  rad-'al  or  "closing”  velocity 
\  of  the  target  with  respect  to  the  missile. 


s 


I 
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•  The  range  of  ”4cpp3er  frequencies”  01  all  the  possible  targets  that  might 
occur  for  this  type  Of  guidance  can  be  determined.  The  band  is  quite  narrow. 

The  velocity  gate  is  swept  through  the  possible  frequencies  until  a  “doppler” 
appears  in  the  gate.  The  antenna  then  "locks  on”  to  the  target  which  has  this 
doppler  frequency  and  obtains  pointing  errors  firm  this  one  target  only, 

“When  the  missile  is  locked  on  to  a  articular  tt-rget,  then  the  direction 
I  of  the  target  Is  known.  Using  the  known  direction  of  the  antenna,  control  in- 

i  formation  for  the  missile  is  derived  so  as  to  produce  a  trajectory  which  will 

j  ■  intercept  the  target. 

\ 

;  The  heming-type  guidance  can  nsqploy  several  alternative  methods  for  intsr- 

;  ception.  The  most  staple  guidance,  of  course,  would  ml newer  the  missile  so 

!  that  it  it  always  headed  toward  the  target.  The  drawbacks  to  this  procedure 

are  immediately  obrlcus.  If  the  target  maneuvers,  the  missile  might  expend  so  1 

.  much  of  its  energy  correcting  its  course,  that  it  is  unable  to  maintain  si  Ti- 

cient  speed  to  reach  Intercept.  Aiao,  if  the  target  nas  a  high  speed,  the 
.  missile  would  have  to  chase  it,  again  erpeadlrgf  a  great  deal  of  energy.  r 


|  Generally,  a  type  of  guidance  is  employed  which,  utilising  the  target 

.  velocity  and  direction,  derive*  a  predicted  intercept  point;  the  "ia-ile  Is 
|  directed  toward  this  paint.  This  method  is  called  proportional  navigation, 
t  Hr  proportional  navigation  we  mean  that  the  rate  cf  change  of  missile  he**  lug 


(that  is,  how  f&e.t  the  missile  must  tars)  must  he  directly  proportional  to  the 
rate  of  rotation  of  the  line-of -sight  fraa  the  missile  to  the  target. 

To  appreciate  the  application  of  this  concept,  consider  the  figure  below. 

In  this  figure,  the  missile  and  target  position  at  e<ja*l  time  Intervals,  t. ,  are 
depicted  in  two  dimensions. 


It  can  he  seen  that  the  line-of-sight  angle,  X  *  does  not  change— hence  the  tar¬ 
get  sad  missile  are  on  a  colli eon  course. 

Ihe  next  figure  shows  the  missile  reaction  if  a  llne-of -eight  rata  exists. 
The  missile  then  attempts  to  correct  its  own  motion  so  that  It  will  follow  a 
constant  heading-angle  that  will  result  in  Intercept. 


Here  we  can  see  that  the  line-of-sight  angle,  X ,  is  changlxg.  The  missile 
must  turn  at  a  rate,  proportional  to  the  llne-of»slght  rate,  in  such  a  mgr  that 
the  change  of  A  is  decreased  to  zero.  When  this  occurs,  the  missile  Is 
i  oriented  along  a  collision  course.  Expressed  mathematically  this  relation  be¬ 
comes: 

»  (On  *  NcOss 
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re  J.e  th-5  rst-ff  of  turn  of  the  missile  heading  (that  is,  the  missile 
velocity  vector)  is  the  rate  of  turn  of  the  llne-of-sight  of  the  missile 

to  the  target,  ard  N  is  the  proportionality  constant.  The  choice  of  N  de¬ 
pends  upon  how  much  response  is  needed  for  the  maneuver  and  whether  the  missile 
can  survive  the  maneuver  without  "breaking.  The  next  figure  indicates  accelera¬ 
tion  as  a  function  of  the  distance  from  the  missile  to  the  point  of  intercept  % 

for  various  values  of  N: 


For  &  medium  speed  target,  N  ■  k  appears  to  be  a  highly  suitable  choice  for 
the  navigation  ratio. 

If  medium  speed  targets  were  the  only  kind  we  expected,  the  choice  of  H  -  4 
for  the  navigation  ratio  would  be  adequate.  Eowever,  we  need  greater  maneuver¬ 
ability  for  high  speed  targets,  but  we  do  not  wish  to  overcorrect  for  slower  tar¬ 
gets.  To  compensate  for  various  target  speeds  an  "effective"  navigational  ratio, 
N',  is  defined  sucl  that 

n'*  vh  ^ 


where  Vc  is  the  closing  velocity  of  the  target  to  the  missile  scd  V*  is  the 
missile  velocity.  Thus,  the  closing  velocity  governs  the  size  of  the  ratio, 

VC/VM,  providing  for  optimaa  maneuverability  for  the  missile. 

The  proportional  navigation  ratio  is  used  to  determine  acceleration  caanands 
to  the  missile.  This  can  be  seen  from  the  following  analysis.  We  knew  that  a 
particle  following  a  curved  path  will  have  a  tangential  velocity, '% ,  equal  to  the 
cross  product  of  the  angular  velocity  and  the  radius -of -curvature  vector.  Thus: 

%=COXR 

Tor  example,  if  the  particle  is  following  a  circular  orbit,  then: 
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The  centripetal,  or  lateral,  acceleration  Is  the  change  of  VT  vith  reepect 
to  the  angular  rate*  Thus,  we  have: 

Al  *  COX  VT  *  Co ft)~  o  o  co  U-6)*R'S* 

•COR  o  o 

or 


A'Vr^S* 

Applying  thla  formula  to  the  total  velocity,  tjj,  of  the  alaalla  and  the 
rate  of  turn  of  the  nlaelle,  ,  we  hertt 

If  we  aatune  that 

V„*V„f  ani  COn‘Q„  F 

<u  -* 

where  I  and  F  are  unit  vector*  for  the  velocity  and  engular  rate  nee- 
peotlvely,  we  oan  write* 
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AL  ~  6l)m  v*  rx £ 


For  the  particular  prvbl^  under  consideration,  we  are  interested  in 
determining  an  acceleration  which  will  change  the  direction  of  the  missile  so 
that  the  line -of -sight  rate  will  go  to  zero.  The  acceleration  desired  will 
he  directed  at  right  angles  to  the  plane  of  the  rotational  vector  and  the  uii  ssil- 
heading  vector.  She  rotational  velocity  (for  the  desired  change)  will  he  the 
rate  of  turn  of  the  missile  heading  OJu  .  The  following  diagram  illustrates 
this: 


For  this  maneuver  (alteration  of  the  heading  angle)  we  can  assume  that  the  angles 
of  attack  are  sufficiently  small  so  that  we  can  replace  the  heading  velocity,  vL, 
ty  the  total  missile  velocity,  vjj.  Our  desired  maneuver,  then  can  he  se*n  from 
the  following  Illustration! 


Calling  A^  the  "command*  acceleration,  we  can  write: 


fa 


Substituting  this  equation  into  the  expression  for  the  proportional 
norlgatlan  ratio  gives: 
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■4^=A/'<y» 

Which  becomes,  using  the  effective  navigational  ratio: 

■$r= 

or 

Thus  from  a  measurement  of  the  sight-line  turning  rate  we  can  derive  lateral 
acceleration  conaands  to  the  missile. 

The  only  acceleration  commands  to  the  missile  are: 

Ak5,A^^m> 

Since,  In  the  previous  derivation,  we  were  only  Interested  In  the  lateral  accele¬ 
ration,  and  since  our  commands  are  only  to  be  given  for  lateral  acceleration,  the 
approximations  we  assumed  are  fully  justified. 


